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SIMULATION  OF  THE  INTERACTION  BETWEEN  AIRDROP 
PLATFORMS  AND  AIRCRAFT  ROLLERS 


INTRODUCTION 

The  C-141  aircraft  is  equipped  with  rollers  on  which  airdrop  platform  loads  ride  during 
loading,  flight,  and  discharge.  The  aircraft  design  places  limits  on  the  load  that  these  rollers 
are  allowed  to  support;  a  single  roller  load  must  not  exceed  1580  pounds*  and  the  total  load 
on  any  row  of  four  rollers  across  the  aircraft  must  not  exceed  2500  pounds.  Before  any 
airdrop  load  can  be  flown  in  die  aircraft  it  must  be  certified  that  these  1g  roller  load  limits 
are  not  exceeded.  This  certification  is  currently  accomplished  by  a  series  of  tests  on  the  roller 
load  test  facility  at  NLABS.  These  tests  require  the  delivery  of  the  airdrop  load  to  the  test 
facility  at  NLABS,  the  rigging  of  the  load  on  the  platform,  and  the  subsequent  testing  which 
requires  20  repetitions  to  obtain  statistical  significance.  Because  of  the  complex  nature  of 
these  tests,  it  was  suggested  that  the  tests  might  be  replaced  with  a  mathematical  simulation. 
Such  a  simulation  would  have  to  be  implemented  on  a  digital  computer  and  would  include 
a  model  of  the  aircraft  roller  system,  a  platform  model,  a  model  of  the  load  to  be  airdropped 
and  the  interaction  of  these  separate  structural  models.  In  this  report  we  describe  some  work 
on  the  development  of  such  a  model.  This  initial  work  is  restricted  to  modeling  the  roller 
system  and  platform  and  their  interaction  under  load.  The  load  used  here  is  a  dead  weight 
load  and  no  attempt  is  made  to  model  real  airdrop  loads.  Two  models  of  the  interaction 
are  treated:  a  deterministic  model  in  which  the  rollers  and  platforms  are  always  in  contact 
and  a  statistical  model  in  which  the  contact  between  the  rollers  and  platform  is  described 
statistically.  In  the  case  of  the  statistical  interaction  model  we  examine  the  influence  of  roller 
stiffness  on  the  distribution  of  load  among  the  rollers.  Earlier  work  on  the  structural  modeling 
of  airdrop  platforms  is  reported  in  reference  1. 

DETERMINISTIC  MODEL 

Initial  considerations  of  the  modeling  of  the  aircraft  roller  system,  the  airdrop  platform 
and  their  interaction  suggested  a  rather  straightforward  linear  structural  model  and  here  we 
describe  this  model,  give  some  computer  results  and  compare  these  results  with  some  test  results. 

Model  Description 

Platform  Model 

The  majority  of  the  structural  model  is  concerned  with  the  air  delivery  platform  and 
here  we  consider  only  the  Type  II  Modular  platform.  This  platform  is  described  in  detail 
along  with  its  rail  system  in  US  Army  NLABS  drawings  11-1-317,  11-1-318,  11-1-319, 


*The  drawings,  nomenclature,  and  computer  programs  describing  the  Type  II  airdrop  platforms 
use  U.S.  Customary  units.  Since  this  report  depends  on  these  data,  these  units  are  used  here 
also.  Table  3  gives  the  conversion  between  U.S.  Customary  and  SI  units. 

'W.  S.  Chang  and  E.  A.  Ripperger;  Stress  and  Deflection  in  Type  II  and  Type  IV  Airdrop 
Platforms;  US  Army  Natick  Labs  Technical  Report  70-56-AD;  Dec  1969.  (AD-711556) 


11— 1—320  and  11— 1—321.  As  illustrated  in  Figure  1 ,  these  platforms  are  aluminum  skin-balsa 
core  sandwich  panels  with  beams  on  the  edges.  The  longitudinal  beams,  those  along  the  short 
sides,  have  provisions  for  the  attachment  of  the  side  rails,  and  these  rails  are  included  in  the 
model.  The  lateral  beams,  those  along  the  long  sides,  serve  as  panel  close-outs.  The  four-foot 
panel  modules  can  be  joined  to  form  various  length  platforms  from  eight  to  twenty-eight  feet 
in  four-foot  increments. 

The  computer  structural  model  is  constructed  using  the  NASTRAN  finite  element  program, 
and  in  Appendix  A  a  FORTRAN  program  to  generate  a  NASTRAN  Bulk  Data  Deck  which 
embodies  this  model  is  presented.  This  FORTRAN  program  includes  documentation  and  model 
data  and  requires  only  the  input  of  an  integer  code  to  specify  the  length  of  the  platform 
for  which  the  model  is  desired.  In  the  following  we  give  a  general  description  of  this  model. 
The  basic  layout  of  the  model  for  the  eight-foot  platform  in  terms  of  the  coordinates,  node 
locations  and  element  types  is  shown  in  Figure  2.  As  shown  in  this  figure,  the  origin  of 
coordinates  is  located  laterally  at  the  center  of  the  platform  and  longitudinally  along  the  first 
row  of  rollers.  This  location  is  common  for  all  platform  lengths,  and  since  the  platforms 
are  located  relative  to  the  rollers  by  the  aircraft  restraint  system,  the  hangover'  on  the  end 
of  the  platform  varies  with  the  platform  length  as  do  the  x  coordinates  of  the  node  locations. 
The  nodal  pattern  is  generally  regular  with  nodes  located  longitudinally  at  roller  positions  and 
laterally  at  roller  positions  plus  three  other  positions  chosen  to  give  reasonable  plate  element 
goemetries.  This  regular  pattern  of  node  placement  is  not  used  at  the  platform  edges  nor 
at  the  junctions  of  panels  where  the  regular  pattern  would  have  resulted  in  unreasonably  narrow 
elements.  In  these  cases  the  nodes  (nodes  36,  37,  ...,  42  for  example)  are  placed  along  the 
edges  of  the  balsa  core-sandwich  panel;  in  fact  the  lines  of  nodes  along  the  panel  edges  define 
the  edge  of  the  balsa  core.  These  lines  of  nodes  do  not  coincide  with  the  centroids  of  the 
beams  but  are  offset.  As  a  result,  this  line  of  nodes  is  also  offset  from  the  rollers,  along 
the  platform  edges  which  contact  the  edge  beams.  The  center  line  of  the  roller  is  outboard 
of  this  line  of  nodes. 

At  the  junctures  of  the  panels  we  have  two  rows  of  closely  spaced  nodes  with  no  intervening 
plate  elements.  For  example  there  are  no  plate  elements  between  the  row  containing  nodes 
36,  37,  ...,  42  and  the  row  containing  nodes  43,  44,  ...,  49.  The  three  inches  between  these 
two  rows  of  nodes  is  the  space  occupied  by  the  lateral  edge  beams  in  the  actual  platform. 
The  two  rows  are  needed  because  the  panels  are  not  joined  with  complete  displacement 
continuity  as  would  be  the  case  if  a  single  row  of  nodes  was  used  along  the  module  juncture. 
The  panels  are  joined  in  the  model  by  equating  the  corresponding  displacements  and  rotations 
at  the  edge  nodes  and  only  the  bending  rotations  at  the  interior  nodes  in  each  of  the  two 
rows.  The  beams  are  modeled  with  the  NASTRAN  beam  element,  CBAR,  which  requires  six 
degrees  of  freedom  at  each  of  its  nodes  so  all  nodes  to  which  beams  are  joining  retain  the 
usual  six  degrees  of  freedom.  The  platform  is  modeled  with  the  plate  bending  element  CQDPLT, 
which  requires  only  three  degrees  of  freedom  per  node,  the  transverse  bending  displacement 
and  two  bending  rotations,  so,  for  all  nodes  to  which  only  plate  elements  are  joined,  only 
these  degrees  of  freedom  are  retained. 

Roller  Model 

The  roller  system  in  the  C-141  aircraft  consists  of  four  rows  of  rollers  laterally  across 
the  cargo  compartment  and  these  rows  are  spaced  every  ten  inches  longitudinally  along  the 
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[~n~|  designates  plaie  element  no.  n 
(n)  designates  beam  element  no.  n 
O  designates  a  node 
^  designates  a  node  connected  to  a  roller 


x.  longitudinal  axis 


Figure  2.  Layout  of  the  element  grid  for  the  platform  model 


aircraft.  The  roller  which  can  be  seen  in  Figure  3  is  a  cylinder  3-3/4  inches  long  with  a 
diameter  of  1-7/8  inches. 

The  roller  is  modeled  as  a  linear  one-dimensional  spring  and  the  actual  NASTRAN  element 
used  is  CELAS2  which  joins  a  specified  node  to  ground  with  a  spring  of  specified  stiffness. 
This  roller  stiffness  was  determined  from  the  test  of  the  roller  setting  on  a  rigid  platform 
as  shown  in  Figure  3.  As  a  result,  the  model  contains  no  influence  of  aircraft  floor  stiffness. 
This  was  done  for  at  least  two  reasons:  data  on  aircraft  floor  stiffness  were  not  readily  available, 
and  in  the  roller  test  facility  the  rollers  rest  on  rather  rigid  supports  and  the  model  will  be 
validated  against  test  data  from  the  roller  test  facility.  The  model  thus  constructed  simulates 
the  roller  test  facility  not  the  aircraft  floor.  Once  the  aircraft  floor  stiffness  is  available,  it 
is  a  simple  matter  to  change  the  model  to  simulate  the  floor.  In  fact,  it  would  be  relatively 
simple  to  have  the  spring  stiffness  nonuniform  with  respect  to  location  on  the  platform. 

Platform-Roller  Interaction 

The  interaction  between  the  platform  and  the  roller  system  is  modeled  by  permanent 
contact  between  the  platform  and  the  rollers.  Since,  in  general,  we  have  a  platform  node 
at  every  roller  location,  this  model  is  accomplished  by  connecting  the  platform  nodes  to  ground 
through  a  spring.  This  is  essentially  enforcing  continuity  between  the  transverse  displacement 
of  the  platform  node  and  the  extension  of  the  spring  connected  to  ground.  There  are  two 
situations  where  the  connection  is  not  so  straightforward:  at  all  the  rollers  along  the  panel 
edges,  and  at  the  juncture  of  the  platform  modules.  In  the  first  case  the  rollers  are  offset 
from  the  nodes  because  the  nodes  are  positioned  to  define  the  edge  of  the  balsa  core  as  discussed 
above.  In  the  second  case  the  nodes  were  offset  from  the  rollers  to  avoid  unreasonably  narrow 
plate  elements.  The  connection  of  the  roller  springs  to  the  platform  nodes  in  the  case  of 
offset  rollers  is  accomplished  by  expressing  the  platform  displacement  over  the  roller  position 
in  terms  of  the  nodal  displacements.  Referring  to  notation  in  Figure  4,  the  continuity  between 
roller  spring  and  platform  is  expressed  as 

ws  =  Wn  +  ?*x  - 

k 

This  can  be  converted  to  a  force  relationship  by  multiplying  by  the  roller  spring  stiffness  K 
to  give 

KWS  =  KWn  +  yK0x  -  xK0y 

The  NASTRAN  code  makes  it  convenient  to  connect  the  node  associated  with  an  offset  roller 
to  ground  with  three  springs  having  stiffnesses  K,  yK,  and  XK  and  associated  with  the  W, 
0X,  and  degrees  of  freedom,  respectively. 

Data  Generation  Program 

The  data  generation  program  presented  as  Appendix  A  is  set  up  to  generate  a  bulk  data 
deck  for  NASTRAN  which  embodies  the  above-described  models.  The  program  requires  input 
of  a  control  parameter  which  indicates  the  platform  length  to  be  modeled  and  some  data 
specifying  the  load  which  is  peculiar  to  the  load  being  modeled  here  and  described  subsequently. 


Test  to  determine  roller  stiffness 


t 

i 
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i  I 

Wn=  DISPLACEMENT  OF  NODE  n 
Ws=  DISPLACEMENT  OF  ROLLER  SPRING 
<#>x>y  =  ROTATIONS  ABOUT  x  AND  y  AXES  AT  NODE  n 
x,  y  =  x  AND  y  OFFSETS  BETWEEN  NODE  AND  ROLLER 


Figure  4.  Nomenclature  for  modeling  the  connection  of  offset  rollers 
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For  any  other  load  distribution  this  part  of  the  data  generation  program  would  need  to  be 
changed.  Since  the  platforms  and  rollers  are  constant,  all  the  data  relative  to  these  models 
is  in  data  contained  in  the  program  and  identified  there.  This  data  includes  the  platform, 
beam  and  roller  stiffness  properties,  and  node  location  properties. 

Model  Performance 

Once  we  have  a  NASTRAN  bulk  data  deck,  it  is  a  relatively  simple  procedure  to  obtain 
a  solution  and  thus  the  roller  load  predictions  of  the  model.  What  is  needed  is  a  NASTRAN 
executive  control  deck  and  case  control  deck  which  are  given  in  Appendix  B  along  with  a 
UNIVAC  EXEC  8  run  stream  to  execute  the  entire  program,  including  generation  of  the  bulk 
data  deck,  on  the  UNIVAC  1106  computer. 

Having  the  model  solution,  it  is  desirable  to  have  a  measure  of  the  performance  of  the 
model  in  predicting  the  roller  loads,  and  the  most  convenient  way  to  get  this  measure  is  to 
compare  the  model  prediction  with  test  results  from  the  roller  test  facility.  To  avoid  the 
complication  of  modeling  an  airdrop  load  such  as  a  truck  or  tank,  a  24-foot  platform  was 
loaded  with  dead  weight  loads  and  put  on  roller  test  facility  to  determine  the  roller  loads. 
The  distribution  of  the  weights  on  the  platform  is  shown  in  Figure  5.  The  nodal  loads  were 
computed  by  assigning  the  platform  area  equally  among  the  nodes  and  the  load  assigned  to 
each  node  is  the  portion  of  the  dead  weight  loads  covering  its  area.  The  data  generation 
program  includes  the  computation  of  the  nodal  loads  based  on  this  procedure  and  the  data 
describing  the  particular  load  distribution  shown  in  Figure  5  are  contained  in  the  file-element, 
ECS* ROLLER. DW— LOADS,  which  is  added  in  the  run  stream  of  Appendix  B.  The  weight 
of  the  platform  was  not  included  in  this  computation,  so  we  should  expect  to  see  a  difference 
between  the  model  results  and  the  test  results  of  a  little  more  than  a  thousand  pounds,  the 
weight  of  the  platform. 

The  results  from  the  24-foot  platform  with  dead  weights  on  the  roller  load  test  facility 
are  presented  in  Table  1  along  with  the  prediction  of  the  roller  loads  given  by  the  model 
described  above.  The  data  in  Table  1  are  the  roller  loads  in  pounds  and  the  table  is  arranged 
so  that  each  row  corresponds  to  a  row  of  rollers  across  the  aircraft;  the  columns  correspond 
to  proceeding  longitudinally  along  the  aircraft.  The  test  results  for  rollers  numbered  30  and 
63  are  blank  because  these  rollers  were  not  operating  or  supporting  the  load  when  the  test 
was  performed.  The  sum  of  the  roller  loads  from  the  test  is  26,446.5  pounds,  and  that  from 
the  model  result  is  25,274.5  pounds  with  the  difference  representing  the  platform  weight  which 
was  not  included  in  the  model.  The  platform  weight  could,  with  relative  ease,  be  added  to 
the  model,  but  here  we  are  interested  more  with  distribution  of  the  roller  loads  than  in 
magnitude,  and  it  is  not  believed  that  this  difference  will  influence  the  conclusion  to  be  drawn 
regarding  the  model  performance  in  predicting  roller  loads. 

We  begin  by  looking  at  the  model  predictions,  and  we  notice  that  they  behave  as  would 
be  expected.  The  load  is  essentially  uniform  across  the  platform  and  each  row  of  rollers  shows 
a  distribution  indicative  of  this  type  of  loading  with  the  load  on  the  two  center  rollers  being 
nearly  equal  and  that  load  being  two  to  three  times  that  on  the  outer  rollers  which  are  also 
nearly  equal.  As  one  looks  along  a  column  of  Table  1  which  corresponds  to  moving 
longitudinally  along  the  aircraft  or  platform,  the  roller  loads  vary  as  they  have  ballast  weights 


Figure  5.  Distribution  of  ballast  weight  on  the  24-foot  platform 


Table  1.  Comparison  of  roller  loads  determined  by  the  model  with  those  from  test. 


a> 

3 

5 


o 

Z 


o^^oi'-.pqiqinoqnoqofl'^oqifi 

iriddddioNffiVirioion 

CO  ^  in  00  ^-CO«-lO*-*-lD«-»- 

co  in  m  mi-  ^  •-  ro  <-  n 


OOlfllfi'JMOl'-N 


r*  co 


O)  N  iO  r*  00 
o>  CO  «-  CO  o> 
lO  CO 


I— 

co 

in 

CO 

CO 

00 

in 

CO 

in 

05 

CO 

CM 

e- 

CO 

CO 

00 

q 

O 

in 

CO 

CO 

cd 

o 

cd 

cd 

05 

in 

IT) 

co 

8 

1“ 

05 

cd 

00 

8 

CN 

l< 

1^ 

cd 

in 

i-L 

CM 

d 

cd 

co 

c» 

co 

CO 

05 

LO 

in 

05 

CO 

00 

CO 

LO 

in 

CM 

r* 

00 

M — 

in 

00 

o 

00 

cm 

*" 

«■“ 

r— 

" 

" 

CM 

00 

00 

a 

00 

o 

o> 

05 

CM 

05 

00 

05 

tn 

o> 

co 

05 

r-> 

05 

00 

O) 

8 

8 

O 

CM 

O 

CO 

o 

s 

in 

o 

8 

o 

8 

8 

o 

- 

CM 

CO 

1 


o 


CM 

CO 

o 

CO 

«- 

05 

T- 

CO 

CO 

e— 

o 

cq 

00 

05 

q 

q 

Mi- 

T“ 

in 

q 

05 

cd 

o 

s 

00 

00 

00 

00 

00 

CD 

t— ^ 

t— ^ 

id 

8 

id 

00 

a> 

O 

Mi¬ 

05 

CM 

CO 

s 

8 

o> 

in 

8 

Pv 

1 — 

T— 

co 

in 

co 

CM 

ce 

CM 

r» 

CM 

in 

CM 

CO 

CM 

0) 

"S 

s 


q 

05 

in 

^r 

o 

o> 

CM 

O) 

o 

05 

CO 

e- 

o 

r** 

CO 

05 

CO 

^r 

00 

CO 

r*. 

CO 

q 

in 

8 

id 

CM 

05 

CO 

r! 

8 

CD 

cd 

CM 

Mf 

ih 

cd 

8 

id 

8 

cd 

o 

id 

CM 

o 

o 

CO 

co 

M* 

co 

in 

CD 

in 

in 

00 

C' 

r*. 

mi- 

o 

00 

CM 

co 

<J> 

CM 

CM 

CD 

*— 

in 

co 

CM 

CO 

*— 

co 

Mf 

CO 

T— 

CO 

in 

«- 

CM 

CO 

CM 

o)0<-CNro^u)(00‘000)0-*Nncifl(or'000)0<-CMP)^ 

incocQioiDiococoiocoiDi^r''r^r-r-'r-r'.r''r'.r-*oooooooQoo 


s 


3 


air'OTrinr>'«fcocM»-ina500>oq'toin^i-'j-oo«-in^- 

idcddoocdcdp^dcdiddricdMrcdbr^odcMCMMrdcdcd 
oooo  *-  oo®T}CMpu5^rco*-ocooocM^-in»-oo»-oo 

CM  T-  CM  CM  *-  MCMfONtM^iniOCOCO^CO 


r-  co 


« 

S 


o 

Z 


o»-CNco^-in<or^ooa>o»-CNCO^-mcor'000>o  —  CMCo^ttn 
nfococjcoooncoron<cv^^'^^«^'^^iDioioioiou) 


2 


8 


cq 

05 

m 

o 

00 

00 

*— 

05 

o 

00 

05 

o 

CO  N  O 

3; 

00 

cd 

Mf 

CM 

cd 

cd 

Mf 

o 

cd 

CT> 

00 

00 

CM 

cd  cd  o 

o> 

CM 

in 

CM 

CM 

CO 

Ml- 

CO 

r— 

in 

05  05 

m 

CM 

CD 

CM 

in 

CM 

?M  r- 

oqooo)N<- 
d  b  cm  cm  cd  id 
o  CM  co 


co'TO^-oocoTrr'OOMrcoooopinoinininooooo>^r«-CM 

gjr^^'cMcdcd^cDoddcdidoddopd'-dcdodddqpoicdcd 

COCOCO.-COOOCO^’-OO'-Tj-'-CD^fO'TCMr^r'CMCMCO'J-CDC'O 


»—  cm  co  ^  in 


ior'000)0>-CNroMrincor'000)pi-CMco^J-in(pr^ooo) 


145.0  37.8  56  488.9  116.2  85  494.1  553.7  114  196.9  58.1 

79.8  0.0  57  254.2  430.5  86  263.5  199.5  115  111.0  6.3 

6.1  2.8  58  42.4  230.3  87  46.2  474.6  116  11.3  0.0 


directly  over  them  or  not.  In  comparing  these  predicted  roller  loads  with  those  obtained  from 
the  test,  we  find  disagreement;  in  fact,  the  two  are  not  in  agreement  in  any  quantitative  way. 
However,  in  a  very  general  way  the  test  results  have  the  same  trends  in  the  distribution  of 
roller  loads  in  that  the  outer  rollers  generally  have  lower  loads  than  the  inner  ones,  and  those 
directly  under  the  ballast  weights  have  higher  toads.  When  the  distribution  of  roller  loads 
is  examined  in  detail,  there  are  some  rather  serious  differences  between  the  actual  test  results 
and  the  expected  distribution.  There  are,  for  example,  numerous  unloaded  rollers,  some  of 
them  being  center  rollers.  In  the  row  of  rolfers  consisting  of  rollers  numbered  7,  36,  65  and 
94  the  left  outer  roller  load  is  652.4  pounds  and  the  left  inner  roller  load  is  3.5  pounds. 
Such  a  combination  of  load  is  completely  contrary  to  that  expected  for  the  type  of  loading 
present  and  in  the  context  of  the  model  described  above.  There  are  other  examples  but 
these  illustrate  the  behavior  found  in  the  test  and  this  behavior  calls  into  question  the 
assumptions  of  the  model,  in  particular,  the  assumption  of  continual  contact  between  platform 
and  rollers  made  in  modeling  the  platform-roller  interaction.  In  the  next  section  we  look 
at  another  assumption  for  this  interaction. 

STATISTICAL  MODEL 

The  results  obtained  with  the  deterministic  model  did  not  agree  well  with  test  results 
and  the  nature  of  the  disagreement  suggests  that  the  model  of  the  interaction  between  platform 
and  the  rollers  was  based  on  a  poor  assumption,  that  assumption  being  that  the  zero  deflection 
positions  of  the  rollers  are  coplanar  and  that  the  platform  is  always  in  contact  with  all  the 
rollers. 

The  data  suggest  that  a  more  accurate  model  would  be  based  on  the  assumption  that 
the  number  of  rollers  in  contact  with  the  platform  is  dependent  on  the  magnitude  of  the 
load.  This  is  caused  by  gaps  between  the  unloaded  platform  and  rollers  due  to  waviness  of 
the  platform  surface  and  unevenness  in  the  roller  heights.  If,  in  the  presence  of  such 
imperfection,  an  unloaded  platform  is  set  on  the  roller  system,  it  will  be  supported  by  only 
a  fraction  of  the  rollers.  As  the  load  is  placed  on  the  platform,  additional  rollers  will  be 
loaded  as  the  platform  deflects.  However,  the  rollers  which  initially  support  the  unloaded 
platform  and  the  loading  sequence  for  the  remaining  rollers  will  be  different  for  each  platform 
and  roller  system  combination.  Since  there  are  many  platforms  and  rollers  systems  in  use, 
it  was  felt  that  a  statistical  model  of  the  platform-roller  interaction  was  needed.  Such  a  model 
complicates  the  analysis  in  two  ways:  first,  a  model  describing  the  relationship  between  the 
rollers  and  the  platform  must  be  developed;  and  second,  the  model  becomes  nonlinear  and 
must  be  solved  incrementally.  In  this  section  of  the  report  we  describe  such  a  model,  its 
solution,  and  present  some  results  given  by  the  model. 

Model  Description 

It  should  be  recognized  that  this  statistical  model  will  be  far  more  complex  than  the 
deterministic  model,  and  because  of  this  increased  complexity  we  used  an  ad  hoc  platform 
structural  model  to  expedite  development  of  a  computational  algorithm  during  the  initial  phase 
of  this  work.  This  ad  hoc  model  has  some  simplification,  is  smaller  in  size  regarding  number 
of  degrees  of  freedom,  it  does  not  have  any  beams  on  the  edges,  and  does  not  make  use 
of  NASTRAN.  We  have,  however,  tried  to  keep  the  general  characteristics  such  as  platform 
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stiffness,  platform  width  and  roller  spacing  the  same  as  the  Type  II  platform-roller  system. 
After  first  describing  this  ad  hoc  model,  we  will  discuss  integration  of  the  NASTRAN  platform 
structural  model  described  above  into  the  statistical  model  described  in  this  section  of  the 
report. 

Platform  Model 

The  addition  of  the  statistical  description  of  the  platform-roller  interaction  makes  the 
problem  nonlinear  and  requires  several  solutions  for  any  loading  case  to  obtain  statistical 
significance.  Both  of  these  factors  greatly  increase  the  amount  of  computation  required;  thus, 
it  seemed  desirable  in  the  development  stages  to  treat  smaller,  less  complex  structural  models 
than  those  described  above  for  the  deterministic  case.  In  addition,  the  nonlinearity  of  the 
problem  made  the  use  of  NASTRAN  unattractive.  As  a  result,  it  was  decided  to  construct 
our  own  ad  hoc  platform  model  using  the  finite  element  described  in  reference  2.  This  is 
a  rectangular  plate  bending  element  with  three  degrees  of  freedom,  transverse  deformation, 
and  the  two  bending  rotations  at  each  of  the  four  nodes.  The  element  has  isotropic  stiffness 
properties  and  does  not  include  transverse  shear  deformation.  In  debugging  the  program  that 
used  this  finite  element,  we  found  that  it  contained  an  error  as  presented  in  reference  2.  This 
error  is  in  matrix  element  (2,  3)  of  submatrix  K, ,  ,  or  matric  element  (8,  3)  of  the  complete 
stiffness  matrix,  and  this  element  should  be  the  negative  of  that  given  on  page  30  of  reference  2. 
An  error  became  apparent  when  we  found  that  the  spectrum  of  the  element  stiffness  contained 
only  two  eigen  values  of  value  zero  when  three,  one  for  each  of  the  rigid  body  degrees  of 
freedom,  are  expected.  Correction  of  the  error  stated  above  resulted  in  the  expected  spectrum. 
The  element  was  also  used  to  solve  the  problem  of  the  simply  supported  rectangular  plate 
loaded  at  its  center  with  a  concentrated  bad.  Satisfactory  results  were  obtained  from  this 
solution  which  was  compared  with  results  from  reference  3.  With  these  checks  on  the  finite 
element,  we  had  confidence  in  it  and  used  it  to  construct  the  platform  model  shown  in  Figure  6. 
Note  that  this  model  has  the  same  geometry  relative  to  platform  width,  node  location,  and 
roller  spacing,  and  the  platform  stiffness  is  the  same  as  that  used  in  the  NASTRAN  plate 
elements.  This  model  differs  from  the  NASTRAN  models  described  above  in  that  it  does 
not  have  any  edge  beams  and  it  is  shorter,  resulting  in  a  reduction  in  the  size  of  the  problem. 

Roller  Model 

The  roller  model  used  here  is  no  different  from  that  used  in  the  deterministic  model. 
It  consists  of  a  simple  scalar  spring  of  stiffness  measured  experimentally  and  given  previously. 
In  the  ad  hoc  model,  all  rollers  were  assumed  to  be  directly  under  nodes,  so  there  was  no 
need  to  treat  offset  rollers. 


2J.  S.  Przemieniecki;  Theory  of  Matrix  Structural  Analysis;  McGraw-Hill  Book  Co.,  New  York; 
1968. 

3S.  Timoshenko  and  S.  Woinowsky-Kreiger;  Theory  of  Plates  and  Shells;  McGraw-Hill  Book 
Co.,  New  York;  1959. 
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Platform  finite  element  model  for  use  in  the  statistical  model 


Platform-Roller  Interaction  and  Solution  Process 

It  is  in  the  platform-roller  interaction  that  the  fundamental  difference  between  the 
deterministic  and  statistical  models  appears.  It  should  be  recalled  that,  in  the  deterministic 
model,  the  platform  and  rollers  were  always  in  contact.  Here,  the  number  of  rollers  in  contact 
with  the  platform  depends  on  the  magnitude  of  the  load.  The  unloaded  platform  is  assumed 
to  sit  on  a  specified  number  of  rollers,  nj,  with  spaces  or  gaps  between  the  remaining  rollers 
and  the  platform.  It  is  in  the  selection  of  the  initial  n;  rollers  and  the  specification  of  the 
size  of  the  gaps  that  the  statistical  nature  of  the  model  lies,  and  these  parameters  are  chosen 
as  follows.  We  obtain  a  set  of  nr  random  numbers  having  a  gamma  distribution  with  specified 
mean  and  standard  deviation.  The  gamma  distribution  is  used  because  it  will  always  give  positive 
gaps.  Each  of  the  random  numbers  is  paired  with  the  integer  name  of  one  of  the  rollers 
resulting  in  an  unordered  array  of  random  numbers  and  a  corresponding  array  of  ordered 
integers.  The  random  numbers  are  then  rearranged  into  ascending  order  with  each  number 
retaining  its  originally  assigned  roller  association.  This  gives  an  ordered  array  of  real  numbers 
and  a  corresponding  array  of  unordered  integers.  Since  the  real  numbers  are  in  ascending 
order,  we  set  the  first  nj  of  them  to  zero,  and  the  rollers  associated  with  these  nj  numbers 
become  the  rollers  on  which  the  unloaded  platform  initially  rests.  The  remaining  nr  —  nj 
random  numbers  are  the  gaps  for  the  rollers  with  which  they  are  associated.  Given  this  model, 
the  solution  for  the  roller  loads  proceeds  incrementally  as  follows:  the  stiffness  matrix  for 
the  platform  is  computed  and  updated  to  include  the  rollers  on  which  the  unloaded  platform 
initially  rests.  This  linear  problem  is  solved  for  the  specified  load.  The  displacements  from 
this  solution  are  compared  with  the  gaps  or  spaces  between  the  platform  and  the  rollers  not 
in  contact  with  the  platform  to  determine  which  roller  will  next  contact  the  platform  and 
what  fraction  of  the  load  is  required  to  cause  contact  to  occur.  Once  this  load  fraction  is 
known,  the  incremental  displacements  associated  with  all  rollers  and  the  incremental  roller  loads 
for  the  rollers  in  contact  with  platform  are  computed  and  saved  as  the  total  displacements 
and  total  roller  loads.  This  completes  the  first  increment  and  we  proceed  to  the  next  by 
updating  the  stiffness  matrix  to  account  for  the  new  roller  that  has  come  into  contact  with 
the  platform.  With  the  solution  of  this  updated  problem,  we  find  a  new  load  fraction;  here, 
however,  we  must  take  into  account  the  deformation  of  the  platform  over  the  rollers  due 
to  previous  increments  as  well  as  the  displacement  of  the  current  increment  and  the  initial 
spaces  between  the  platform  and  rollers.  Given  the  new  load  fraction,  the  incremental  roller 
loads  and  displacements  are  computed  and  used  to  update  the  total  roller  loads  and 
displacements.  After  each  increment,  a  check  must  be  made  to  see  if  any  rollers  have  been 
unloaded,  since  the  rollers  must  be  in  compression  or  unloaded.  If  any  roller  is  unloaded, 
a  new  load  fraction  which  causes  unloading  to  occur  is  computed  and  the  total  roller  loads 
and  displacements  are  again  updated  using  this  load  fraction.  The  stiffness  matrix  is  then 
updated  either  by  removing  a  roller,  if  one  has  become  unloaded,  or  by  adding  a  roller  for 
the  next  one  to  be  contacted  by  the  platform.  This  process  then  continues  until  either  all 
the  load  is  acting,  at  which  point  we  have  the  total  roller  loads  for  the  rollers  in  contact 
with  the  platform  and  perhaps  some  unloaded  rollers,  or  until  all  rollers  are  in  contact  with 
the  platform,  in  which  case  the  final  increment  is  assigned  the  load  fraction  required  to  apply 
all  the  remaining  load.  In  either  case,  the  result  is  the  roller  load  distribution  associated  with 
the  specified  statistical  model  of  the  platform-roller  interaction.  Such  a  distribution  is  computed 
several  times,  as  specified,  and  the  statistics  of  the  individual  roller  loads  and  the  sum  of  the 
roller  loads  in  each  row  of  four  are  computed  and  printed  as  output.  A  FORTRAN  program 
to  carry  out  this  computational  process  is  presented  in  Appendix  C. 


Evaluation  of  the  Model 


The  success  of  a  statistical  model  such  as  the  one  described  above  depends  on  one's 
knowledge  of  the  statistic  of  the  gaps  between  the  rollers  and  the  platform.  At  present  such 
knowledge  is  totally  absent.  Not  only  do  we  not  know  the  statistical  parameters  of  the  gamma 
distribution  used,  but  we  do  not  even  know  if  the  gamma  distribution  is  suitable  for  the 
representation  of  the  phenomenon.  Neither  of  these  deficiencies  represents  a  problem  from 
a  model-building  or  program  point  of  view,  since  both  the  statistical  model  and/or  the  statistical 
parameters  can  be  changed  with  ease  and  with  no  impact  on  the  rest  of  the  solution  process. 
However,  from  the  point  of  view  of  model  evaluation,  these  deficiencies  are  crippling  and  no 
real  evaluation  can  be  made.  The  elimination  of  these  deficiencies  would  require  an  extensive 
experimental  program  to  measure  the  roller-platform  gaps.  Such  an  experimental  program  should 
include  variations  in  aircraft  and  platforms  in  order  to  get  meaningful  data. 

Regarding  such  data,  observation  of  the  platform  used  in  the  roller-load  tests  described 
here  indicated  a  waviness  in  the  platform  that  might  best  be  described  in  terms  of  spectral 
quantities  such  as  period  and  amplitude,  with  these  spectral  parameters  being  specified  in  a 
statistical  fashion.  Additionally,  it  appeared  that  the  platform  waviness  was  the  most  significant 
part  of  the  imperfect  matching  in  comparison  to  the  unevenness  of  the  rollers  in  the  roller 
test  facility.  The  unevenness  of  the  rollers  in  an  actual  C— 141  aircraft  may  present  quite 
a  different  picture.  Until  such  data  are  obtained,  the  statistical  roller  load  model  described 
here  is  incomplete,  and  little  in  the  way  of  evaluation  can  be  done. 

In  constructing  this  model,  it  became  apparent  that  the  magnitude  of  the  roller  stiffness 
was  so  large  that,  with  the  allowed  loads,  the  rollers  would  undergo  displacements  on  only 
one-  or  two-tenths  of  an  inch.  It  is  easy  to  believe  that  unevenness  in  the  rollers  or  waviness 
in  the  platform  of  this  order  of  magnitude  exists,  and  thus  that  the  occurrence  of  unloaded 
rollers  is  not  only  possible  but  likely.  This  suggests  that  making  the  rollers  much  more  flexible 
would  eliminate  the  occurrence  of  unloaded  rollers  and  thus  improve  the  distribution  of  load 
among  the  rollers.  To  demonstrate  this  we  used  the  statistical  model  described  above  to  run 
a  series  of  calculations  in  which  we  varied  the  roller  stiffness  and  looked  at  two  measures 
of  the  load  distribution,  the  maximum  load,  and  the  maximum  sum  of  the  load  on  the  four 
rollers  in  a  row.  In  these  calculations  we  used  a  load  of  10,000  pounds  divided  equally  among 
all  the  nodes  in  the  model  and  the  statistical  model  used  had  a  mean  and  standard  deviation 
of  0.002;  when  the  mean  and  standard  deviation  are  equal,  the  gamma  distribution  reduces 
to  the  exponential  distribution.  The  results  of  these  calculations  are  shown  in  Table  2  where 
we  present  the  mean  and  standard  deviation  of  the  maximum  roller  load  and  the  maximum 
row-sum  of  roller  loads.  In  addition,  Table  2  also  contains  the  mean  and  standard  deviation 
of  the  differential  maximum  roller  load  and  the  differential  maximum  row-sum  of  roller  loads. 
These  differential  quantities  are  computed  by  obtaining  the  difference  between  these  parameters 
given  by  the  statistical  model  and  by  a  direct  solution  with  all  rollers  in  contact  and  the  full 
load  on  the  platform.  Thus,  the  differential  quantities  are  measures  of  departure  from  what 
might  be  called  "the  ideal  situation."  The  statistics  of  these  roller  load  parameters  are  based 
on  twenty  determinations  of  the  roller  load  distribution  with  the  statistical  model.  The  results 
in  Table  2  show  that  whenever  statistical  behavior  is  included,  the  maximum  roller  loads  are 
greater  than  those  occurring  in  the  ideal  case,  but  that  the  variation  from  the  ideal  decreases 
as  the  roller  stiffness  decreases.  In  fact,  with  the  data  in  Table  2,  it  cannot  be  claimed  in 
a  statistical  sense  that  the  differential  maximum  row-sum  of  roller  loads  is  nonzero  when  the 
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Table  2.  Variation  of  Maximum  Roller  Loads  With 

Roller  Stiffness 

Maximum  Roller  Load,  lb 

Total 

Differential 

Roller 

Stiffness,  Ib/in 

Mean 

Std.  Dev. 

Mean 

Std.  Dev. 

r  ■ 

5,000 

-400.622 

1.86 

-5.486 

1.86 

*  '■ 

?  *• 

25,000 

451.240 

10.91 

-26.489 

10.91 

fc’ 

50,000 

479.083 

14.41 

49.693 

14.41 

75,000 

-506.005 

26.92 

-72.56 

26.92 

£ 

100,000 

-519.157 

20.72 

-83.05 

20.72 

>7: 

125,000 

-557.971 

43.86 

-119.92 

43.86 

150,000 

-553.227 

54.84 

-113.69 

54.84 

Maximum  Row-Sum  of  Roller  Loads,  lb 

i  * 

Tflrai 

Differential 

• 

Roller 

t-. 

Stiffness,  Ib/in 

Mean 

Std.  Dev. 

Mean 

Std.  Dev. 

<* 

5,000 

-1453.10 

3.47 

-2.49 

3.47 

1 

25,000 

-1467.29 

18.31 

-27.94 

18.31 

50,000 

-1495.60 

27.89 

-61.28 

27.89 

75,000 

-1521.64 

36.24 

-88.11 

36.24 

100,000 

-1525.77 

34.39 

-90.05 

34.39 

125,000 

-1555.27 

52.46 

-117.80 

52.46 

1 

i 

■ 

150,000 

-1545.86 

37.03 

-106.95 

27.03 
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roller  stiffness  is  5,000.  This  is  the  result  of  the  estimate  of  the  mean  being  less  than  twice 
the  standard  deviation.  It  is,  in  fact,  less  than  the  standard  deviation.  Thus  we  see  that 
as  the  roller  stiffness  decreases,  the  roller  load  distribution  approaches  the  ideal.  The  maximum 
roller  stiffness  for  which  data  are  presented  in  Table  2  is  approximately  the  stiffness  of  the 
rollers  now  in  use,  and  this  causes  a  25%  increase  in  the  maximum  load  for  the  rather  small 
imperfection  of  a  0.002  inch  average  gap  between  rollers  and  unloaded  platform.  It  seems 
clear  that  more  flexible  rollers  would  be  an  improvement  relative  to  the  roller  load  distribution 
problem,  since  this  decreases  the  magnitude  of  the  maximum  roller  loads  and  makes  the 
distribution  of  roller  loads  nearer  the  ideal  distribution. 

It  seems  appropriate  to  say  a  few  words  about  the  amount  of  computing  required  to 
get  solutions  from  this  statistical  model.  We  begin  by  noting  that  this  ad  hoc  model  has  49 
nodes  with  three  degrees  of  freedom  at  each  node  for  a  total  of  147  degrees  of  freedom. 
Thus  for  each  increment  in  the  solution  a  system  of  equations  of  order  147  must  be  solved. 
Of  these  49  nodes  in  the  model,  28  are  associated  with  rollers.  Of  these  28  we  have  assumed 
in  the  calculation  presented  that  eight  are  initially  in  contact  with  the  platform,  leaving  20 
to  be  added  incrementally,  one  on  each  increment.  Typically,  two  or  three  additional  increments 
are  required  because  of  unloading  of  rollers,  so  to  get  one  sample  set  of  roller  loads,  22 
increments  are  required.  To  obtain  significant  measures  of  the  statistical  behavior  of  the  roller 
loads,  we  have  been  generating  twenty  samples,  and  these  calculations  take  85  minutes  of  central 
processor  time.  Thus,  the  amount  of  computer  time  even  for  this  rather  small  model  is  rather 
large.  The  computing  time  can  be  reduced  by  generally  cleaning  up  the  computation  process 
and  also  by  carrying  out  a  condensation  of  the  platform  structure  before  the  incremental  process 
is  stated.  This  condensation  would  remove  all  the  bending  rotations  and  the  transverse 
displacements  not  associated  with  rollers.  For  this  ad  hoc  model  such  a  condensation  would 
reduce  the  problem  of  order  147  to  a  problem  of  order  28  and  this  would  have  a  vast  impact 
on  the  computation  time.  Now,  what  is  of  more  interest  is  what  this  might  tell  us  relative 
to  computation  for  a  statistical  model  with  the  NASTRAN  platform  models.  These  NASTRAN 
models  have  respectively  40,  56,  76,  96,  116,  and  132  rollers  on  the  8,  12,  16,  20,  24,  and 
28-foot  platforms,  so  that  with  condensations  the  order  of  the  problems  to  be  solved  would 
be  equal  to  the  number  of  rollers  plus  one-half  that  many  again  to  account  for  rotations  needed 
for  offset  roller  connections.  It  then  appears  that  the  ad  hoc  model  with  147  degrees  of 
freedom  is  near  the  higher  limit  of  the  NASTRAN  models  relative  to  the  order  of  the  problem 
to  be  solved.  The  other  factor  that  affects  computation  time  is  the  number  of  iterations 
required,  and  this  depends  on  the  number  of  rollers  and  the  number  of  rollers  initially  in 
contact  with  the  platform.  Since  the  ad  hoc  model  has  28  rollers  with  eight  in  contact  initially, 
all  the  NASTRAN  platform  models  will  require  more  iterations  than  the  ad  hoc  model.  From 
these  considerations,  it  appears  that  with  the  NASTRAN  models  the  computing  time  might 
be  significantly  shorter  with  the  small  platforms  and  significantly  larger  with  the  large  platforms 
than  was  experienced  with  the  ad  hoc  model.  This  is  all  based  on  the  assumption  that  the 
statistical  model  in  which  a  small  percent  of  the  rollers  are  initially  in  contact  with  the  platform 
is  correct  and,  as  we  have  stated  above,  we  have  not  been  able  to  evaluate  this  model.  If 
it  should  be  shown  that  a  statistical  model  is  one  in  which  a  large  percent  of  the  rollers  are 
in  contact  initially,  this  would  greatly  reduce  the  computation  time  by  reducing  the  number 
of  iterations  required  for  solution. 

Integration  of  the  NASTRAN  platform  models  into  the  statistical  solution  is  discussed 
in  Appendix  D. 


CONCLUDING  REMARKS 


In  a  study  of  the  mathematical  modeling  of  the  loads  on  C141  aircraft  cargo  rollers  resulting 
from  airdrop  loads,  two  models  were  constructed:  a  deterministic  model,  and  a  statistical 
model.  The  deterministic  model  proved  to  be  inadequate  because  it  failed  to  represent  the 
imperfection  in  the  platform-roller  interface.  It  is  believed  that  this  imperfection  has  a  statistical 
character.  To  include  this  character,  a  statistical  model  was  developed  and  tested  successfully 
in  a  mathematical  sense.  However,  because  of  lack  of  data,  it  was  not  possible  to  specify 
the  statistical  character  of  the  interface,  and  thus  the  model  could  not  be  evaluated  in  the 
physical  sense.  Further  data  are  needed  to  determine  both  the  qualitative  and  the  quantitative 
character  of  the  platform-roller  interface.  The  statistical  model  was  used  to  show  that  the 
use  of  more  flexible  rollers  would  improve  the  roller  load  distribution  by  decreasing  the 
significance  of  the  imperfection  relative  to  the  deflections  of  the  rollers. 
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Table  3.  Conversion  table.  US  Customary  units  to  SI  units. 


To  convert  US  Customary 


Quantity 

US  Customary 

SI  Units 

units  to  SI  units 
multiply  by 

Mass 

pounds  mass 

kilograms 

0.455 

Force 

pounds  force 

newtons 

4.45 

Length 

inch 

meter 

0.0254 

foot 

meter 

0.305 

yard 

meter 

0.91 

Area 

square  inch 

square  meters 

6.45  x  10'4 

square  foot 

square  meters 

0.093 

Volume 

cubic  inches 

cubic  meters 

1.64  x  10  s 

cubic  feet 

cubic  meters 

0.0283 

Density 

pounds  per 
cubic  inch 

kilograms  per 
cubic  meter 

2.77  x  104 

ounces 

grams  per 
square  meter 

34 

Tension 

pounds  per  inch 

newtons  per  meter 

176 

Moment  of 
Inertia 

(inches)4 

(meters)4 

4.1  x  10'7 

Modulus  of 
Elasticity 
and  Stress 

pounds  per 
square  inch 

newtons  per 
square  meter 

6.89  x  103 
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APPENDIX  A 

FORTRAN  Program  to  Generate  NASTRAN 
Bulk  Data  Deck 
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fKECSDlNO  PAOX  BUMUMOC  FILMED 


DATA-GEN 


1.  COMP I LER(DlAG*3) 

2.  OIMENSION  IXS<6) ,IXF(6) ,XT< 164) ,YT(7> 

3.  OIMENSION  LLC ( 6 ) , KKC ( 6 ) . I ROL (53 ) , YROL ( 18) 

4.  DIMENSION  XT! (55) ,XT2( 109) 

5.  DIMENSION  ZL0FSL(3) ,ZL0FSR(3) .ZTOFSBO) ,ZT0FST(3) 

6.  DIMENSION  QPLPP(5)  ,QBPRL< '•) ,QBPRT(5) .QMAT1 (4) ,QMAT2(4) 

7.  EQUIVALENCE  ( XT  1 ( 1 ) , XT ( 1 ) ) , ( XT2 ( 1 ) , XT ( 56 ) ) 

8.  C 

9.  c  definition  of  the  DATA  arrays 

10.  c  yt-  y  coordinates  of  platform  GRID  POINTS 

11.  C  XT1+XT2=XT-  X  COORDINATES  of  platform  grid  points 

12.  C  IXS.IXF-  THE  START  AND  FINISH  LOCATIONS  IN  XT  FOR 

13.  C  the  different  LENGTH  platforms 

14.  C  LLC-  STARTING  POINTS  IN  IRCL  FOR  DATA  FOR  DIFFERENT 

15.  C  LENGTH  PLATFORMS 

16.  C  IROL-  POSITIVE  VALUES  DENOTE  POSITION  NUMBERS  ALONG 

17.  C  LENGTH  OF  PLATFORM  AT  WHICH  NODES  ARE  NOT  CONNECTED 

18.  C  TO  ROLLER  SPRINGS  NEGATIVE  VALUES  DENOTE  POSITIONS 

19.  C  WHERE  CONNECTION  TO  ROLLER  IS  THRU  AN  OFFSET,  REF  BY  LLC 

20.  C  YROL-  X  OFFSETS  BETWEEN  GRID  POINTS  AND  ROLLER  SPRINGS 

21.  C  REF.  BY  KKC  AND  NEG. VALUES  OF  IROL 

22.  C  KKC-  STARTING  POINTS  IN  YROL  FOR  DIFFERENT  LENGTH  PLATFORMS 

23.  C  ZLOFSOFFSETSFOR  LONGITUDINAL  BEAM  BENDING  AXIS 

24.  C  FROM  NOOAL  AXIS 

25.  C  ZTOFS-  OFFSET  FOR  TRANSVERSE  BEAM  BENDING  AXIS 

26.  C  FRON  NODAL  AXIS 

27.  C  QPLPR-  SECTION  PROPERTIES  FOR  PLATFORM  PLATE 

20.  C  QBPRL-  SECTION  PROPERTIES  FOR  LONGITUDINAL  BEAM 

29.  C  QBPRT-  SECTION  PROPERTIES  FOR  TRANSVERSE  BEAM 

30.  C  ROLSTF  «  XOFS-  ROLLER  STIFFNESS  AND  OFFSET  AT 

31 .  C  EDGE  OF  PLATFORM 

32.  C  QMAT1,  QMAT2-  MATERIAL  PROPERTIES  FOR  AL  $  BALSA  WOOO 


33. 

DATA 

IXS/1 ■ 13,32.56,86 

122/ 

34. 

Data 

IXF/12,31 ,55,85,121 , 164/ 

35. 

data 

YT/-50 . 53 1 .-32.8 

-14.8.0. 

.14.8.32 

B, 

36. 

50.531/ 

37. 

DATA 

XT  1 /  -3.75  ,  0. 

.  10.0 

,  20.0 

,  30.0 

,  41.125, 

38. 

44.125,  50.0 

.  60.0 

,  70.0 

,  BO.O 

,  89.0  , 

39. 

-2.75  ,  0.0 

.  10.0 

,  20.0 

,  30.0 

,40.0  , 

40. 

42.125,  45.125.  50.0 

.  60.0 

,  70.0 

,  80.0  , 

41  . 

90.0  ,  93.0 

.100.0 

.110.0 

.120.0 

,130.0  , 

42. 

138.875, 

43. 

-6.75  .  0.0 

.  10.0 

,  20.0 

.  30.0 

,  38.125, 

44. 

41.125.  50.0 

,  60.0 

,  70.0 

,  80.0 

.86.0  , 

45. 

89.0  ,100.0 

,110.0 

,120.0 

,130.0 

,133.875, 

46. 

136.875,140.0 

.150.0 

.160.0 

,170.0 

,181.75  / 

47. 

data 

XT2/  -0.75  ,  10.0 

.  20.0 

,  30.0 

.  40.0 

,  44.125, 

48. 

47.125.  50.0 

,  60.0 

,  70.0 

,  80.0 

,92.0  , 

49. 

95.0  ,100.0 

.110.0 

.120.0 

,130.0 

,139.875, 

50. 

142 . B75 , 150.0 

,160.0 

,170.0 

,180.0 

,137.75  , 

51  . 

190.75  ,200.0 

.210.0 

,220.0 

.230.0 

,235.625, 

52. 

-4.75  ,  0.0 

.  10.0 

,  20.0 

,  30.0 

,  40.125, 

53. 

43.125,  50.0 

,  00.0 

,  70.0 

,  80.0 

,  88.0  , 

54. 

91.0  ,100.0 

,  1  ••  0 . 0 

,120.0 

,130.0 

, 135.075, 

55. 

138.875,150.0 

,150.0 

.170.0 

.180.0 

,183.75  , 

56. 

186.75  ,190.0 

,200.0 

,210.0 

.220.0 

,231 .625, 

57. 

234.625,240.0 

.250.0 

,260.0 

,270.0 

.279.5  . 
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DATA-GEN 


58. 

-3.75  ,  0.0  ,10.0 

,  20.0  .  30.0  ,  41.125, 

59. 

44.125,  50.0  ,  60.0 

,  70.0  ,  80.0  ,  89.0  , 

60. 

92.0  ,100.0  .110.0 

,120.0  ,130.0  ,136.875, 

61 . 

139.875,150.0  ,160.0 

,170.0  ,100.0  ,184.75  , 

62. 

187.75  ,190.0  ,200.0 

,210.0  ,220.0  ,230.0  , 

63. 

232.625,235.625.240.0 

,250.0  ,260.0  ,270.0  , 

64. 

280.5  ,283.5  ,290.0 

,300.0  .310.0  ,320.0  , 

65. 

328.375/ 

66. 

DATA 

LLC/ 1,5,10,18,28, 40/ 

67. 

DATA 

KKC/1 ,1,3,6,10,15/ 

68. 

OATA 

IROL/  1.  -6,  7,-12,  1. 

7,  8,  14.  19,  1,  6, 

69. 

-7.  12,-13,  IB,  19,-24,  -1,  6,  7,-12,  13,-18 

70. 

19,  24,-25,  30.  1,  - 

6,  7,  12,-13,  18,-19,  24 

71  . 

25.-30,  31,-36,  1,  - 

6,  7.-12,  13,  18,-19,  24 

72. 

25,  31.  32,-37,  38,  43/ 

73. 

DATA 

YROL/-I .125,  1.0  ,-1.125, 

1.0  ,-1.75  ,  0.75  . 

74. 

-2.0  ,  0.125.-0.75  , 

-0.125,-1.0  ,  1.125, 

75. 

-1.625,  0.5  ,-1.125, 

1.0  ,  0.125,  0.5  / 

76. 

data 

ZLOFSL , ZLOFSR/O . 0 , -1 .795.-0 

.2648,0.0,1 .795,-0.2648/ 

77. 

data 

ZT0FSB.ZT0FST/-1 .0525,0.0,0 

.0,1 .0525,0.0,0.0/ 

78. 

DATA 

QPLPR/O. 163.2.5,0.0248, 1 .3, 

-1.3/ 

79. 

DATA 

QBPRL/4 .46 ,2.69,5.93,8.62,0 

.0/ 

80. 

DATA 

<JB°RT/0. 6687,0. 148,0.6683,0 

.8163,0.0/ 

81  . 

data 

ROLSTF , XOFS/3000 .0,0. 269/ 

82.  DATA  QMAT1/10000000. 0,4000000. 0,0. 33, 0.098/ 

83.  DATA  QMAT2/1 .0,24900.0,0. 1 .0.006/ 


84. 

85. 

86. 

87. 

88. 

89. 

90. 
91  . 

92. 

93. 

94. 

95. 

96. 

97. 

98. 

99. 
100. 
101  . 
102. 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110. 
Ill  . 
112. 

113. 

114. 


C 

C  THE  PARAMETER  'LENGTH'  SPECIFIES  PLATFORM  LENGTH  CONSIDERED 

C 

C  LENGTH  123456 

C  PLATFORM  LENGTH  8  12  16  20  24  28 

C 
C 

C  GENERATION  OF  GRID  CARDS  CONTAIN  NODE  NUMBER 

C  ANO  NODAL  COORDINATES  IN  THE  IMPLIED  GLOBAL 

C  COORDINATE  SYSTEM 

C 

C  JYN=  THE  NUMBER  OF  NODES  IN  THE  Y  DIRECTION,  ACROSS  THE 

C  WIDTH  OF  THE  PLATFORM 

C 

10  FORMAT ( ) 

WRITE(13,1012) 

1012  FORMAT (  1  BEGIN  BULK') 

RE AD( 5 ,10)  LENGTH 
IS=IXS( LENGTH) 

I F* I XF( LENGTH) 

JYN.7 
I NODE-O 
Z*0.0 
JW=C 

DO  100  I* IS,  IF 
X*XT ( 1 ) 

UW  =  ,!W*1 

UWT.JV 

UBEAMcO 

IF<(I.EO.IS).OR.(I.EQ.IF))  JBEAM.1 
I F( (LENGTH. ED. 6) . AND. < JW.GT.26) )  JWT-JWT-1 
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DATA-GEN 


115. 

1 16. 

117. 

118. 

119. 

120. 
121  . 
122. 

123. 

124. 

1 25 . 

126. 

127. 

128. 

129. 

130. 
131  . 

132. 

133. 

134. 

135. 

136. 

137. 

138. 

139. 

140. 
141  . 

142. 

143. 

144. 

145. 

146. 

147. 

148. 

149. 

150. 
151  . 

152. 

153. 

154. 

155. 

156. 

157. 

158. 

159. 

160. 
161  . 
162. 

163. 

164. 

165. 

166. 

167. 

168. 

169. 

170. 
171  . 


128 

100 

1000 

1014 

C 

C 

C 

C 

C 

C 

c 

c 


1008 


C 

C 

C 

c 


IF( (LENGTH. EQ.2) .AND. (dW.GT.2) )  dWT.dWT-1 
JTEST«MOO( JWT.6) 

IF(( JTEST . EO.O) .OR. (JTEST .EQ.1 ) )dBEAM*1 

00  100  d*1,dYN 

L6EAM.0 

IF( ( J.EQ.1 ) . OR . ( J . EQ. JYN ) )  LBEAM* 1 
Y=VT(d) 

1  NOPE* I  NODE  +  1 

IF( (OflEAM.EO. 1 >. OR. (LBEAM. EQ.1 ) )  GO  TO  128 
NR  I TE ( 1 3 , 1 000 )  INODE, X.Y.Z 
GO  TO  100 

WRI TE ( 1 3 , 1 01 4 )  INODE, X.Y.Z 
CONTINUE 

FORMAT ( ■  GRID1 ,3X. 18 . 8X. 3FB . 3 , 1 3X . '126' ) 

FORMAT ( 1  GRID' ,3X, IB.8X,3F8.3) 

GENERATION  OF  DUMMY  GRID  POINT  CARDS  FOR  USE  IN 
DEFINING  BEAM  LOCAL  COOORDINATE  SYSTEMS 

1GR I D 1  FOR  LONGITUDINAL  BEAMS  ON  LEFT 
IGRI02  FOR  TRANSVERS  BEAMS  ON  BOTTOM 
IGRI03  FOR  LONGITUDINAL  BEAMS  ON  RIGHT 
IGRI04  FOR  TRANSVERSE  BEAMS  ON  TOP 

IGRI D1 • INODE  +  1 
IGRI02* INODE  +2 
J  GRI 03= I  NODE  *  3 
I GR 1 04* I  NODE  +  4 
X*350 . 0 
Y=  1  (»0 . 0 

WRITE( 13, 1008)IGRID1 , X.Y.Z 

FORMAT (  ‘  GRI0'3X.  I8.3X.3FS..-  1  OX .' 1 23456 '  ) 

X=-50.0 

WR!TE( 13, 1008) IGRI 02, X.Y.Z 
Y=-1 00 . 0 

NR  1 TE ( 1 3 , 1 008 ) I GR 103, X.Y.Z 
X  =  350.0 

WRITE! 13. 1008) IGRI 04, X.Y.Z 

GENERATION  OF  COOPT  CARDS  FOR  THE  QUAORI LATERAL  BENDING 
elements  for  the  platform  model 


IEIM*0 

lQP*lXFi LENGTH)  -  IXS(LENGTH) 

JYN1  *v)YN  -  1 
DO  101  1*1 . IQP 

DO  101  d«1 . dYNI 

Iw*  i 

1F( ( LENGTH. SO- 6) . AND. ( I . GT . 25 ) ) I W* tN- 1 
IF< ( LENGTH. 00-2) . AND. ( I .GT .1 ))IW=IW-1 
I5K!P  =  V,0P(  IW.6) 

I r ( ISKIP.EO.O)  GO  TO  101 

d 1 =  d  +  (  1-1  ) *dYN 

02  =  d 1  +  ' 

d3=d 1 +JYN 

d4  =  J3+ 1 

IELM*IF.LM+1 

WHITE! 13.1001 )  IELM.d1.d2.d4.d3 
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OATA-GEN 


172. 

173. 

174. 

175. 

176. 

177. 

178. 

179. 

180. 
181 . 
182. 

183. 

184. 

185. 

186. 

187. 

188. 
189. 
1*0. 
191  . 

192. 

193. 

194. 

195. 

196. 

197. 

198. 

199. 

200. 
201  . 
202. 

203. 

204. 

205. 

206. 

207. 

208. 

209. 

210. 
211  . 
212. 

213. 

214. 

215. 

216. 

217. 

218. 

219. 

220. 
221  . 
222. 

223. 

224. 

225. 

226. 

227. 

228. 


101  continue 

1001  FORMAT!'  CODPLT  * , I8.7X, *  1 ' ,4I3,5X, '0.0' ) 

C 

C  GENERATION  OF  PROPERTY  CAROS  FOR  CQDPLT  PLATE  ELEMENT 

C 

WRITE! 13. 1002)  OPLPR 

1002  FORMAT ('  PQDPLT  ' ,  7X,  •  1  ’  .  7X,  • 1 ' , F8.3. 71, *2 ' ,4FB. 3) 

C 

C  GENERATION  OF  MATERIAL  PROPERTY  CARD 

C  QMAT1  6  0MAT2  CONTAIN  YOUNG'S  MODULUS . SHEAR 

C  MOOULUS.  POtSSON'S  RATIO  AND  MASS  DENSITY 

C 

C  MATERIAL  PROPERTY  for  BENDING. REF.  ON  PQDPLT 

C 

WRI TE{ 1 3 , 1 003 )QMAT 1 

1003  FORMAT ( 1  MAT  1  * . 7X , • 1 ' . 2E8 . 2 . 2F8 . 3 ) 

C  MATERIAL  PROPERTY  FOR  TRANSVERSE  SHEAR,  REF. ON  PODPLT 

C 

WRITE! 13,1 004 )  QMAT2 

1004  FORMAT ( 1  MAT  1  ' , 7X , 1 2 ‘ , 2E8 . 2 , 2F8 . 3 ) 

C 

C  GENERATION  OF  THE  CBA3  CARDS  FOR  THE  BEAMS  ALONG  THE 

C  PLATFORM  EDGES 

C 

c  LONGITUDINAL  BEAMS 

C 

1 CONT  *  1 00 

iqpi.iqp  ♦  i 

DO  112  1*1 ,IQP 

JL1-JYN*( 1-1 )  +  1 
JL2*JL1  +  JYN 
JRIrJLI  ♦  JYN  -  1 
JR2*JR 1  ♦  JYN 
1 PBARL* 1 
I E  LM  *  I E  LM  +  1 
I CONT* I CONT  ♦  1 

WRITE! 13, 1009)IELM,IPBARL.dL1 ,JL2, IGR1D1 , I CONT . I CONT , ZLOFSL , ZLOFSL 
I  ELM* I  ELM  +  1 
I CONT* I CONT  +  1 

WRITE! 13, 1009  HELM, IPBARL , JR2 , JR1 , IGRI03, ICONT, ICONT .ZL0FSR.2L0FSR 
1009  FORMAT! '  CBAR ' . 3X , 518, 23X , 1 20' , J4/ , J4 , 1 9X , 6F8 . 3) 

112  CONTINUE 
C 

C  TRANSVERSE  BEAMS 

C 

LCNP2=LENGTH  +  2 
I PBART=2 

DO  113  I.1.LENP2 
INC*6 

MOEG* ( 1-1 )«INC»dYN 

IF!! LENGTH. EG. 6) .AND. ( I.GT. 5) )  MBEG=M8EG  +  JYN 

IF( 1  LENGTH. EG. 2) .AND. ( I.GT. 1 ))  MBE3*MBEG  *  JYN 

00  1 13  J*1 , JYN1 

MBEG*M3EG  +  1 

JB1 *MBEG 

JS2* JB1  +  1 

JT1 * J81  -JYN 
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DAIA-GEN  ***** 


229. 

JT2*«JT1  +  1 

230. 

1 F( I . EQ. LENP2 )  GO  TO  114 

231  . 

1  ELM* I ELM  ♦  1 

232. 

1COVT.-ICONT  +  1 

233. 

WR I TE( 1 3, 1009)1  ELM, IP3ART.JB1 . JB2. IGRI02. ICONT , ICONT .ZTOFSB.ZTOFSB 

234. 

114 

CONI  IN'JE 

235. 

1F(I.EQ. 1 )  GO  TO  115 

236. 

1  ELM* I ELM  ♦  1 

237. 

ICONT* ICONT  +  1 

238. 

WRITE(13,1009)IELM,IPBAST,JT2,JT1 , IGRID4. ICONT, ICONT, ZTOFST.ZTOFST 

239. 

115 

CONTINUE 

240. 

1 1 3 

CONTINUE 

241  . 

C 

242. 

C 

GENERATION  OF  THE  PBAR  CAROS  SPECIFYING  THE  BEAM 

243. 

C 

SECTION  PROPRTIES,  QBPRL  6  QBPRT  CONTAIN  AREA, 

244. 

C 

2  AREA  MONENTS.  TORSIONAL  INERTIA  ANO  NONSTRUCTURAL  MASS 

245. 

C 

246. 

C 

247. 

c 

FOR  LONGITUDINAL  BEAMS 

248. 

c 

249. 

I D*  1 

250. 

MID*1 

251  . 

WRITE! 13. 101 0)10, MID, QBPRL 

252. 

1 

1010 

FORMAT! '  PBAR1 ,3X,2I8,5F8.3) 

253. 

c 

254. 

c 

FOR  TRANSVERSE  BEAMS 

255. 

c 

256. 

10*2 

257. 

WRITE! 13. 101 0)10. MID, QBPRT 

258. 

c 

GENERATION  OF  SCALAR  SPRING  ELEMENTS  REPRESENTING  THE 

259. 

c 

ROLLERS,  THESE  HERE  ARE  CONNECTING  THE  PLATFORM  TO  GROUND 

260. 

c 

rolstf*roller  stiffness 

261  . 

c 

XOFS*X  OFFSET  OF  EDGE  NODES  FROM  ROLLER  CENTER  LINE 

262. 

c 

263. 

120*0 

264. 

I  F ! x  2Q • EO . 0 ) GO  TO  778 

265. 

KC=KKC( LENGTH) 

266. 

LC*  LLC( LENGTH) 

267. 

DO  106  1*1 , 10P1 

268. 

YY  *0.0 

269. 

INODE*! 1-1 ) • JYN  -1 

270. 

IF(I.NE.1R0L(LC))  go  TO  102 

271  . 

LC*LC  +  1 

272. 

GO  TO  106 

273. 

102 

CONTINUE 

274. 

IF(I.EQ.IABS(!ROL(LC)))  GO  TO  104 

275. 

GO  TO  105 

276. 

104 

YY*ABSI YROL(KC) ) 

277. 

LC*  LC  +  1 

278. 

K  C  *  ►  C  *  1 

279. 

105 

CONT INUE 

260. 

00  109  J«1.UYN,2 

281  . 

AX *0.0 

282. 

I F ( J.EO. 1 )  XX*XOFS 

283. 

I'!.. EO. JYN)  XX.XOFS 

284. 

inode* inode  ♦  2 

285. 

ZSW*  ROLST  F 
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286. 

287. 

288. 
289. 

290  • 

291  . 

292. 

293. 

294. 

295. 

296. 

297. 

298. 

299. 

300. 
301  . 

302. 

303. 

304. 

305 . 

306. 

307. 

308. 

309. 

310. 
311  . 

312. 

313. 

314. 

315. 

316. 

317. 

318. 

319. 

320. 
321  . 

322. 

323. 

324. 

325. 

326. 

327. 

328. 

329. 

330. 
331  . 

332. 

333. 

334. 

335. 

336. 

337. 
238. 

339. 

340. 
341  . 
342. 


zkphx«rolstf*yv 
zkphy*rolstf*xx 
I  ELM* I ELM  ♦  1 
ID0F*3 

WRITE( 13. 1005)  I  ELM, ZKW, INODE, IDOF 
IF(YY.EO.O.O)  GO  TO  103 
JELM*1ELM  ♦  1 
ID0F*5 

WRITE! 13. 1005  HELM. ZKPHX, INODE, I  DOF 
103  CONTINUE 

IF(xx.EO.O.O)  GO  TO  109 
I  ELM* I  ELM  ♦  1 
I  DO  F  =  4 

WRITE! 13, 1005)  I  ELM, ZKPHY, INODE. IDOF 

109  CONTINUE 

106  CONTINUE 

1005  FORMAT! ‘  CELAS2  ' , 18 , E8. 3 . 218 ) 

778  CONTINUE 
C 

C  SINGLE  POINT  CONSTRAINTS  TO  FIX  THE  PLATFORM 

C  IN  ITS  PLANE 

C 

WR ITE(13,101S) 

1015  FORMAT!1  SPC1 ' ,10X, '1 1 ,6X, '12* ,7X, 1 1 ’ ,7X, ' 7 1 ) 

C 

C  GENERATE  MULTI-POINT  CONSTRAINT  CARDS  THAT  SIMULATE  THE 

C  JUNCTURE  BETWEEN  PLATFORM  SECTIONS 

C  BENDING  ROTATIONS  CONTINUOUS  AT  ALL  NODES  ALONG  JOINT 

C  TRANSVERSE  DISPLACEMENT  CONTINUOUS  AT  PLATFORM  EDGES 

c 

I  SET* 1 

DO  107  1=1 .LENGTH 

INC  =  6 

NBEG«I*INC«JVN 

1 F ( ( LENGTH . EO . 6) . AND. ( I .GT .4 ) )NBEG*NBEG  ♦  JYN 
I F ( ( LENGTH . EQ ■ 2 ) . AND. ( I . GE . 1 ) )NBEG*NBEG  ♦  JYN 
DO  107  J= 1 , JYN 
NBECi=N8EG  +  1 
N3EC.PJ  =  NBEG  ~  JYN 

IF((J.E0.1 ) .OR. ( J.EQ.JYN) )  GO  TO  110 
GO  TO  111 

110  IDOF  =  3 

WRITE!  13. 1007) I  SET, NBEG, I  DCF .NBEGPJ , IDOF 

111  I D0F*4 

WRITE! 13, 1007) I  SET, NBEG, I  DCF. NBEGPJ, IDOF 
I  DO  F  *  5 

107  WRITE! 13, 1007) I  SET .NBEG, I  DOF, NBEGPJ, IDOF 
1007  FORMAT!1  MPC ' . 4X , 3 18, 5X , ' 1 . 0 ' , 2 IB ,4X , ' -1 .0 ' ) 

C 

C  GENERATION  OF  LOADING  CARDS 

C  NODAL  LOADS  FOR  DEAD  WEIGHT  LOADING  TEST  ON  THE 

C  24  FOOT  PLATFORM  IN  ROLLER  TEST  FACILITY 

C 

F  T  =  0 . 0 

READ! 5 .10)  NLOAD 
IF(NLOAD.EQ.O)  GO  TO  129 
DO  117  I  LOAD*  1, NLOAD 


r 

[ 

- 

: 

OATA-GEN  •  **»• 

1 

343. 

WR I T  E  T  6 , 10)  FTS 

344. 

RE ADIS. I0)0LQA0,X1 ,X2,Y1 ,Y2 

345. 

FTS. 0.0 

346. 

I X* 1 XS( LENGTH )  -  1 

) 

£ 

347. 

DO  1 18  1«1 . I0P1 

I 

r 

348. 

NBCGM  1-1  )«JVN 

349. 

1  1  X  -  I  X  ♦  i 

r 

350. 

IF ( 1 .EQ. 1 )  GO  TO  119 

351  . 

GX 1 =0 . 5’ ( XT ( I IX )  ♦  XT ( I IX- 1 )  I 

352. 

GO  10  120 

353. 

119  GXI.XM  IIX)  -  1.5 

354. 

120  IF( ! .EQ. I QP1 )  GO  TO  121 

,  ■ 

355. 

CX2r0.5*(XT( IIX)  ♦  XT ( 1 IX+ 1 ) ) 

356. 

GC  10  122 

357. 

121  GX2.XT  1  1  IX)  ♦  1.5 

358. 

122  CONTINUE 

359. 

I F ( GX2 . LT . XI )  GO  TO  118 

360. 

1 F ( (GX2.GE. XI ) .ANO. 1GX2.LE.X2) . AND. (GX1 . LE.  XI ) )  XL*GX2-X1 

361  . 

IF( (GX2 -GE.X2) . AND. (GX) . LE.X1 ) )  XL.X2-X1 

362. 

IF( (GX2. LE.X2) .ANO. (GX1 .GE.X1 ) )  Xl*GX2  -  GX1 

363. 

IF( (GX2.GE.X2) .ANO. (GX1 .GE.X1 ) . AND . ( GX 1 . LE . X2 ) )  XL-X2-GX1 

364. 

I F ( GX 1 . GT . X2  )  GO  TO  117 

’ 

365. 

00  123  J.I.dVN 

366. 

N0EG=N13EG  ♦  1 

367. 

1 F 1 J . EQ . 1 )  GO  TO  124 

. 

368. 

GY  1 sO . 5* (YT(J)  ♦  YT(J-I)) 

369. 

GO  70  125 

370. 

124  GY1.-52.281 

371  . 

125  IF(,t.EQ..1YN)  GO  T0126 

372. 

&Y2=0.5*( YT(d)  1  VT ( J+1 ) ) 

373. 

GO  TO  127 

374. 

126  GY7-.52.2B1 

375. 

127  CONTINUE 

376. 

IF(GY2. LT . Y1 )  GO  TO  123 

377. 

IF( ( GY2 . GE . Y 1 ) . AND. (GY2.LT. Y2) . AND . ( GY  1 . LT . Y 1 ) )  YL-GY2-Y1 

378. 

I F ( (GY2. LT. Y2) . ANO. (GYt.GT.YI))  YL-GY2-GY1 

379. 

IF( (GY2.GT.Y2) .ANO. ( GV1 .GT . Y 1 ) . AND . ( GY  1 . LT . V2 ) )  Y L« Y2-GY1 

380. 

IF(GY* .GT.V2)  GO  TO  118 

381  . 

FOPCE=XL*YL»QLOAD 

| 

382. 

ISET.t 

383. 

icoro.o 

384. 

VI *0.0 

385. 

V7*0 . 0 

1  1 

386. 

V3* 1  .0 

387. 

WR 11  El  13.101 1 ) I  SET ,N8EG, I  CORO, FORCE, VI ,V2,V3 

388. 

1011  FORMAT  1  FORCE  ' ,3IB,4F8.3j 

389.. 

FT. FT  .  FORCE 

390. 

■  FTS.FTS  +  FORCE 

391  . 

123  CONTINUE 

392. 

lie  continue 

393. 

117  CONI INUE 

394. 

WRITE(6,1C)  ft 

395. 

129  *7:tE(  -.3,10)3) 

396. 

1013  FORMAT) 'ENDDATA' ) 

397. 

STOr 

398. 

END 
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APPENDIX  B 

Procedure  of  Execution  of  Deterministic  Model 
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APPENDIX  B 


Procedure  of  Execution  of  Deterministic  Model 

•  s 

Presented  here  is  the  UNIVAC  Exec  VIII  runstream  for  execution  of  the  data  generation  - 
program  which  yields  the  BULK  DATA  DECK  in  TPF$.DATA  and  the  subsequent  execution 
of  NASTRAN.  This  is  followed  by  the  executive  and  case  control  deck  for  NASTRAN  which 
is  called  by  the  runstream  statement  @ADD,P  R.  EX— CACONT.  The  final  item  is  the  data 
used  by  the  data  generation  program  to  put  the  loads  in  the  bulk  data  deck. 

RUNSTREAM 

•ASQ-AX  ECS*  ROLLER 
@USE  R.(  ECS*  ROLLER 

@ASG,T  13. 

@FOR,S  R.DATA— GEN.DG 
@XQT 
1 

@ADD,P  R.DW— LOADS 
@ELT,I  TPFS.DAT A 

@ADD,P  13. 

@FREE  13. 

@QUAL  NASTRAN 

@ASG,AX  NASTRAN*NASTRAN 

@HDG,N 

@XQT  NASTRAN*  NASTRAN.  LI  NK1 

NASTRAN  HICORE  =  85000  CONFIG  *  14  SYSTEM(34)  *  2 
@ADD,P  R. EX— CACONT 
@ADD,P  TPF$.DATA 

@ADD,P  NASTRAN*NASTRAN.NASTRAN 
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EXECUTIVE  AND  CASE  CONTROL  DECKS 

ID  ROLLER, CHECK 

APP  DISP 

SOL  1,0 

TIME  60 

DIAG  22 

CEND 

TITLE  =  ROLLER  LOAD  SIMULATION 

SUBSTITLE=  8  FT  PLATFORM,  DEAD  WEIGHT  LOADS 

ELFORCE  =  ALL 

STRESS  =  ALL 

OLOAD  =  ALL 

DISP  =  ALL 

ECHO  -  BOTH 

LOAD  =  1 

MPC  =  1 

SPC  =  1 
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CMDGD'JOXJI'fcCjK) 


«»•■>  DW-LOADS 


1  .  13 

1.850,7.05.19.05,-46.681,1.319 
1.826,19.05,31 .05,-46.081,1.919 
1 .832, 6.65,29.65, 1 .919,49.919 
1 .862,54.82,78.82,-48.081 ,-0.081 
1 .8 04,53.02,77.02,-0.62,47.38 
1 .858,99.20,123.20,-45.62,2.38 
1.804,98.30,122.30,1.42,49.42 
1 .849, 148.88,172.88,-45.62,2.39 
1.812,148.78,172.78,-1  .382,46.618 

11 .  1.804, 192.95,216.95,-44.610,3.38 

12.  1.875, 196.25,220.25,2,72,50.718 

13.  1 . 804. 243. 72, 267. 72, -45. 41B, 2. 582 

14.  1 .812,244.02,268.02,1.12,49. 12 
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APPENDIX  C 


Program  for  the  Solution  of  the  Statistical  Model 

The  FORTAN  computer  program  for  solution  of  the  statistical  model  consists  of  a  main 
program  called  MAIN2  and  this  program  uses  six  subroutines  peculiar  to  it.  Subroutine  PLATE 
generates  the  platform  finite  element  structural  model  and  stores  the  stiffness  matrix  which 
is  the  embodiment  of  this  model  on  file  3  in  banded  symmetric  format.  In  carrying  out 
this  computation,  the  subroutine  PFE3  is  called  from  subroutine  PLATE  to  generate  the 
12  x  12  element  stiffness  matrices.  The  subroutine  WMTX  is  conditionally  called  by  subroutine 
PLATE  to  print  as  output  the  element  stiffness  matrices  generated  by  subroutine  PFE3.  These 
matrices  are  printed  if  the  control  parameter  (PRINT  has  the  value  zero.  Subroutine  STGAP2 
is  called  by  the  main  program  to  generate  a  set  of  random  numbers  having  a  gamma  distribution 
that  are  used  to  represent  the  spaces  or  gaps  between  the  platform  and  rollers.  These  random 
numbers  are  paired  with  an  integer  which  specifies  the  roller  with  which  it  will  be  associated. 
Once  paired  with  these  integers,  the  random  numbers  are  placed  in  ascending  order  while 
retaining  the  originally  paired  integer.  This  ordering  is  accomplished  by  a  call  of  subroutine 
AORDER  from  STGAP2.  The  load  vector  is  generated  by  subroutine  PLOAD  which  is  called 
from  the  main  program.  The  computation  in  PLOAD  uniformly  distributes  the  total  load 
specified  as  input  over  all  the  nodes,  both  those  associated  with  rollers  and  those  not.  In 
addition,  the  program  uses  three  general  purpose  subroutines  from  the  IMSL;  GGAMT  to 
generate  the  random  numbers  with  a  gamma  distribution  is  called  from  STGAP2,  LEQ1PB 
to  solve  the  linear  equations  in  band  symmetric  form,  and  BECORI  to  compute  the  statistics 
of  the  roller  loads  after  all  the  incremental  solutions  have  been  completed.  The  subroutine 
GGAMT  requires  a  starting  seed  number  for  which  we  use  the  time  obtained  by  a  call  to 
the  UNIVAC  function  ERTRAN. 

The  input  to  the  program  in  the  order  read  is  as  follows: 

A  —  platform  finite  element  length 

B(1),B(2),B(3)  -  platform  finite  element  widths 

D  -  platform  bending  stiffness 

XNU  —  Poisson’s  Ratio  of  platform  material 

P  -  total  magnitude  of  the  uniform  load  on  the  platform 

ROLK  -  roller  stiffness 

I  PR  I  NT  -  print  control  parameter  used  in  subroutine  PLATE 

I  PR  I  NT  =  0  Print  element  stiffness  matrix 
IPRINT  #  0  Do  not  print 

NEL  —  number  of  elements  in  platform  finite  element  model 
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NUMRO  —  number  of  rollers  initially  in  contact  with  platform 


NT  —  number  of  solutions  for  statistical  accuracy 

NR  —  number  of  rollers 

GM,GSD  —  mean  and  standard  deviation  of  the  Gamma  distribution 
associated  with  the  roller-platform  gaps 

All  of  this  input  is  read  in  free  field  format  and  an  example  of  the  input  is  presented  as 
part  of  the  runstream  for  execution  of  the  program. 

There  are  two  sets  of  output  for  this  program;  one  of  which  is  an  abbreviated  set  and 
is  sent  to  logical  unit  6,  the  standard  FORTRAN  print  file.  This  output  contains  only  a 
brief  summary  for  each  increment  in  the  calculation  containing  the  roller  numbers  in  contact 
with  the  platform,  the  next  roller  to  make  contact,  the  load  fraction,  and  the  sums  of  the 
roller  loads  and  the  incremental  roller  loads.  Following  this  output  is  a  listing  of  the  means 
and  standard  deviations  of  each  of  the  roller  loads  and  of  the  maximum  roller  load  and  maximum 
sum  of  rows  of  four  roller  loads.  The  other  set  of  output  contains  details  from  each  incremental 
solution  including  in  addition  to  that  in  the  summary  output;  a  listing  of  rollers  and  their 
associated  gaps;  incremental  solution  vector;  incremental  roller  displacements  and  loads;  and 
total  roller  displacements  and  load.  This  detail  output  is  sent  to  logical  unit  7  which  is  a 
permanently  assigned  file.  If  the  summary  output  reveals  some  difficulty  in  the  computation 
process  the  detailed  printout  is  then  available  for  use  in  debugging. 

A  listing  of  the  FORTAN  source  programs  follows,  and  a  UNIVAC  Exec  VIII  runstream 
for  execution  of  the  program  and  a  sample  data  set  are  also  presented. 


! 
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IMA  I  N2  »*»«  = 


1  . 
2. 
3. 
4  . 

5. 

6. 

7. 

8. 
9. 

10. 
1 1  . 
12. 

13. 

14. 

15. 

16. 
17. 
16. 

19. 

20. 
21  . 
22. 
23. 
24  . 

25. 

26. 

27. 

28. 

29. 

30. 
31  . 

32. 

33. 

31. 

35. 

36. 

37. 

38. 

39. 

40. 
41  . 

42. 

43. 
44  . 

45. 

46. 

47. 
48- 

49. 

50. 
51  . 

52. 

53. 

54. 

55. 

56. 

57. 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c. 

c 

c 


c 

c 

c 

c 


CCMPIl£R(DIAG*3) 

DIMENSION  NODE(50) ,GAP(50) 

DIMENSION  I RONO ( 28 ) , B ( 3 ) , GK ( 1 47 , 27 ) , BL ( 1 47 , 1 ) 

DIMENSION  X ( 28 ) , F ( 28 ) , XINC ( 28 ) , F INC( 28 ) 

DIMENSION  XT (28) ,FT(28) ,R LA (50, 30) 

DATA  IRONO/1 .7,13,19,22.28.34,40.43,49,55. 

61,64,70,76,82,65.91,97,103,106, 
112,118,124,127.133.139, 145/ 

READ  DATA  FOR  PLATE  GEOMETRY  AND  STIFFNESS 
ASB(K)  ELEMENT  LENGTHS  AND  WIDTHS 
O&XNU  PLATE  STIFFNESS  AND  POISSON'S  RATIO 
NEL  NUM8ER  OF  ELEMENTS 
P  LOAD  MAGNITUDE 
NUMRO  NUMBER  OF  ROLLERS 
RQLK  ROLLER  STIFFNESS  CONSTANT 
IRONO  ARRAY  DEFINING  THE  RELATION  BETWEEN 
ROLLER  NO.  AND  ThEIR  DOF  NUMBER 
I RCNO ( I  )  =  THE  DOF  NO. FOR  ThE  ROLLER  I 

READ(5,10;A,B,D,XNU.P,ROLK 
RE AD (5, 10) I  PR  I  NT. NEL. NUMRO 
READ ( 5 . 1 0 ) NT , NR , GM , G3D 
WRITE(6,23)A,8,0. XNU.NSL.ROlK.P 
VkR !  T  E  (  7 . 23  )  A  ,  3 ,  D .  XNU  .  NE  l  .  ROLK .  P 
23  rORKAT( 1H1 ,//2X. 'ELEMENT  LENGTH' ,F10. 5/ 

2X.  *  ELEMENT  AiD'HS'  .3F10.5/ 

2  X .  'STIFFNESS'  .C15.5/ 

2X,  ’ POISSONS  RAT JC'  , FI  0.5/ 

2X , 1  NO .  OF  PLATE  ELEMENTS'  , 15/ 

2X, 'ROLLER  STIFFNESS' ,E15. 5/ 

2  X ,  ' LOAD'  ,£15.5) 

GENERATE  UNCONSTRAINED  GLOBAL  STIFFNESS  MATRIX 

CALL  PLAT-E(A.B.D,XNU,N£L.  I  PRINT) 

10  FORMAT  |  ) 

NT  -  NUMBER  OF  TESTS  DR  REPI1 ITIONNS  TO  BE  MADE 
TO  GST  A I N  STATISTICAL  SIGNIFICANCE 
NR=NUM6ER  QP  ROLLERS  SJPpGRTING  PLATFORM 
GMsKEAN  OF  ThE  RO lLE R-p LA TFURM  GAPS 
G5D- ST aNDARC  DEVIATION  OF  THE  ROLLER-PLATFORM  GAPS 
NUMRO=NO.  CP  ROLLERS  INITIALLY  IN  CONTACT 
with  the  platform 


WRI TE( 6 , _ J  >nR , GM . GSD , N r , Nt MRO 
WRITE(7.20)NR, CM . 3SD . NT , NUMPO 
20  FORMAT)/, '2X, 'NO.  OF  ROLLERS', 15/ 

?X, 'ROLLER  GAP  MEAN',E15.5/ 

2X , ' ROLLER  GAP  STD . DE V . ' , E 1 5 . 5/ 

2 X  ,  'NO.  OF  TESTS'  . 15/ 

2X , "NO .  OF  ROLLERS  INITIALLY  IN  CONTACT', 15) 

BEGIN  A  LOOP  ON  THE  NUMBER  OF  REPETITIONS  OR 
TESTS.  GET  NCDEE  NUMBERS  FOR  ROLLERS  ON  WHICH 
THE  PLATFORM  INITIALLY  SETS 
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MA1N2  *■*»= 


58. 

NNT  *0 

59. 

DO  100  IR»1 .NT 

60. 

WRITE(6, 26)18 

61  . 

WRITE(7, 26)18 

62. 

26 

FORMAT! 1H1 ,//2X, ‘BEGINNING  OF  TEST  OR  REPETITION ‘ , 14! 

63. 

NNTxNNT  +  1 

64. 

C 

65. 

C 

NNT  RELATES  TO  THE  NUMBER  OF  TESTS  OR 

66. 

C 

REPETITIONS-  WHEN  STIFFNESS  MATRIX  IS 

67. 

C 

FOUND  SINGULAR  NNT  IS  DECREASED  SO 

68. 

c 

SO  IT  IS  THE  TRUE  NO.  OF  TESTS 

69. 

c 

70. 

KEY  =  0 

71  . 

CALL  STGAP2(N0DE,GAP,NR,GM,GSD  ) 

72. 

DO  135  NZ  =  1  . NUMRO 

73. 

135 

GAP ( NZ ) =0 . 0 

74. 

WRITE(6. 19) 

75. 

WR 1 TE ( 7 , 19) 

76. 

19 

F0RMAT(//2X, 'ROLLER  NOS.  &  GAPS') 

77. 

WRITE (6, 1 2) (NODE! IP) .GAP! IP) , IP=1  ,NR) 

78. 

WRI TE ( 7 , 12) (NODE! IP) .GAP! IP) , I P= 1 , NR ) 

79. 

1  2 

FORMAT (2X.I5.E12.5) 

80. 

c 

81  . 

c 

INITIALIZE  PARAMETERS  FOR  INCREMENTAL  SOLUTION 

82. 

c 

83. 

INCNUM=0 

84. 

INCNO=0 

85. 

S.cRAC  =  0 . 0 

86. 

DO  1 10  K- 1 ,NR 

87. 

X ( K ) =0 . 0 

88. 

110 

F ( K ) =0 . 0 

89. 

NRC=NUMRO 

90. 

109 

CALL  PL0AD(tL,P,49) 

91  . 

I NCNUM= I NCNUM  +  1 

92. 

I NCNO  = I NCNO  +  1 

93. 

WRITEf  6,27)1 NCNUM 

94. 

WRITE (7, 27) I NCNUM 

95. 

27 

FORMAT! 1H1 ,//2X, ‘BEGINNING  OF  INCREMENT  * , 14 ) 

96. 

REWIND  3 

97. 

READO, 2001 )((GK(IP.dP),IP«1, 147), JP=1 ,27) 

98. 

2001 

FORMAT (2X.6E20.12) 

99. 

DO  101  1  =  1  ,  NRC 

100. 

L=NODE( I ) 

101  . 

L  L  =  I RONO ( L ) 

102. 

101 

GK( LL,27)=GK< LL.27)  +  ROLK 

103. 

C 

104. 

C 

OBTAIN  SOLUTION 

105. 

C 

106. 

CALL  LE01PB(GK,147,26, 147.BL.147.1 .5.21 .22, I ER) 

107. 

WR I TE ( 7 , 1 0 ) I ER 

108. 

WRITE (6. 33) 

109. 

WRI *E; 7. 33) 

110. 

33 

FCF».AT(//2X,  'ROLLERS  IN  CONTACT  FOR  THIS  INCREMENT') 

Ill  . 

WR1TE(6,34)(N0DE!IP),IP*1 .NRC) 

112. 

WRITE(7,34)( NODE ( I P ) , I P* 1 ,NRC) 

113. 

34 

FORMAT ( 815 ) 

114. 

IF(IER.EO.O)  GO  TO  lit 
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mm* 

MA1N2 

aati 

■  3 

172. 

FT(IX)=0.0 

173. 

XT( IX)=0.0 

174. 

XINC(IX)cO.O 

175. 

120 

FINC ( IX ) =0 .0 

176. 

KC=0 

177. 

RL=0 . 0 

178. 

RLINC=0 . 0 

179. 

FMAX=0 . 0 

180. 

DO  105  J«1,NR 

181  . 

L=N0DE ( J ) 

182. 

LL= I R0N0 ( L ) 

183. 

XINC(L)  =  FRAC*BL(LL,1  ) 

184. 

XT ( L ) =X ( L )  +  X INC ( l. ) 

185. 

IF(J.GE.NRC)  GO  TO  106 

186. 

FINC(l)=-FRAC*BL(LL, 1 )*ROLK 

187. 

FT ( L ) =F ( L )  +  FINC(L) 

188. 

RL  =  RL  +  FT(  L) 

189. 

RLINC=RLINC  +  FINC(L) 

190. 

C 

191  . 

C 

CHECK  TO  SEE  IF  INCREMENT  HAS  CAUSED 

192. 

C 

ROLLER  UNLOADING,  IF  UNIOAOIN  OCCURS  COMPUTE 

193. 

C 

SMALLEST  LOAD  FRACTION  TO  CAUSE  ZERO 

194. 

c 

LOAD  ON  AN  UNLOADED  ROLLER 

195. 

c 

196. 

IF(FT(L).LE.FMAX)  GO  TO  106 

197. 

UFRAC=-F( L)/(BL(LL. 1 )*R0LK) 

198. 

KC=KC  +  1 

199. 

IF ( KC. EO. 1 )  GO  TO  122 

200. 

IF(UFRAC.GE.UTEM)G0  TO  106 

201  . 

122 

UTEM*UFRAC 

202. 

JUNL=J 

203. 

LUNL=L 

204. 

KEY  =  1 

205. 

106 

CONTINUE 

206. 

105 

CONTINUE 

207. 

UFRAC=UTEM 

208. 

c 

209. 

c 

ADJUSTMENT  OF  MOOEL  FOR  ROLLER  UNLOADING 

210. 

c 

211  . 

IF(KEY.EO.O)  GO  TO  121 

212. 

I F ( INCNO.GT. 1 )GO  TO  124 

213. 

c 

214. 

c 

ROLLER  UNLOADED  ON  FIRST  INCREMENT,  INTERCHANGE 

215. 

c 

unloaded  roller,  an  initial  ROLLER,  WITH  NEXT 

216. 

c 

ROLLER  TO  BE  CONTACTED  AND  REPEAT  FIRST  INCREMENT 

217. 

c 

218. 

INCNO=INCNO  -  1 

219. 

NODE( JSAV)=LUNL 

220. 

NODE ( JUN  L ) *  I SAV 

221  . 

NRC  =  NRC“ 1 

222. 

SFRAC=SFRAC-FRAC 

223. 

WRI TE( 6, 35 )  LUNL.ISAV 

224. 

WRIT  E ( 7 , 35 )  LUNL.ISAV 

225. 

35 

FORMAT(//2X, ‘ROLLER  UNLOADED  ON  FIRST  INCREMENT1/ 

226. 

2X, 1  INTERCHANGE  ROLLERS ', 14 , '  AND1, 14) 

227. 

WRITE(7, 19) 

228. 

WRITE(7,12)(NODE(IP),QAP(IP),IP-1,NR) 
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115. 

NNT.NNT  -  1 

1 1B. 

GO  TO  100 

117. 

111 

CONTINUE 

118. 

WRITE(7,14) 

119. 

14 

FORMAT(///2X, 'SOLUTION  VECTOR') 

120. 

00  115  IW  >1 .49 

121  . 

IWW=( IW-1  )»3 

122. 

115 

WRITE(7,13)IW. BL( I WW+1 , 1 ) . BL( 1WW+2 . 1 ) • BL( IWM+3, 1 ) 

123. 

13 

FORMAT (2X, I4.3E15.6) 

124. 

C 

125. 

C 

COMPUTE  LOAD  FRACTION  FOR  CURRENT  INCREMENT 

126. 

C 

127. 

NRC=NRC+1 

128'. 

IF(NRC.GT.NR)GO  TO  116 

129. 

TEMs 100.0 

130. 

00  102  I=NRC,NR 

131  . 

L=NODE( I ) 

132. 

L  L=  I RONO ( L ) 

133. 

IF(L.LT.O)  GO  TO  104 

134. 

I F ( BL( LL. 1 ) . LE.O .0 )G0  TO  104 

135. 

FR AC= ( GAP ( I)-X(L))/BL(LL.1) 

136. 

IF(I.EO.NRC)  GO  TO  103 

137. 

IF(FRAC.GE.TEM)  GO  TO  104 

138. 

103 

T  EM=  FRAC 

139. 

I S AV=  L 

140. 

JSAV=I 

141  . 

104 

CONTINUE 

142. 

102 

CONTINUE 

143. 

frac=tem 

144. 

SFRTsSFRAC 

145. 

SFRAC^SFRAC  +  FRAC 

146. 

1 F ( SFRAC . LE . 1 .0)  GO  TO  107 

147. 

FRAC= 1.0-  SFRT 

148. 

SFRAC=1 .0 

149. 

GO  T0  107 

150. 

116 

FRAC=1 .O-SFRAC 

151  . 

SFRAC=1 .0 

152. 

WRITE(5.28)FRAC,SFRAC 

153. 

WRI1E(7,2B)FRAC,SFRAC 

154. 

2e 

FORMAT ( //2X , 'FINAL  INCREMENT'/ 

155. 

2X, ' INCREMENTAL  LOAD  FRACTION' ,E15. 5/ 

156. 

2X.  'TOTAL  LOAD  FR ACT  ION '  , El  5 . 5 ) 

157-. 

107 

CONTINUE 

158. 

SFRT=SFRAC 

159. 

WR I T  E ( 6 , 22 )  ISAV.FRAC.SFRAC 

160. 

WRITE(7,22)  ISAV.FRAC.SFRAC 

161  . 

22 

FORMAT (//2X, 'NEXT  NODE  TO  CONTACT', 15/ 

162. 

2X, ' INCREMENTAL  LOAD  FRACTION' ,E15. 6/ 

163. 

2X,  'TOTAL  LOAD  FRACTION ', El  5 . 6) 

164. 

IF(hEY.EQ.O)  GO  TO  129 

165. 

I F ( KEY . EQ . 2 )  NODE( JUNL)»-NODE( JUNL) 

166. 

KEY-0 

167. 

129 

CONTINUE 

168. 

C 

169. 

C 

update  deflection  and  roller  loads 

170. 

C 

171  . 

DO  120  IX-1 , NR 
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229. 

230. 
231  . 

232. 

233. 

234. 

235. 

236. 

237. 

238. 

239. 

240. 
241  . 

242. 

243. 

244. 

245. 

246. 

247. 

248. 

249. 

250. 
251  . 

252. 

253. 

254. 

255. 

256. 

257. 

258. 

259. 

260. 
261  . 
262. 

263. 

264. 

265. 

266. 

267. 

268. 

269. 

270. 
271  . 

272. 

273. 

274. 

275. 

276. 

277. 

278. 

279. 

280. 
281  . 
282. 

283. 

284. 

285. 


KEY«0 
GO  TO  109 
124  CONTINUE 

I F ( LUNL . NE . LROLA )  GO  TO  128 
NODE ( JUNL ) *-NOOE( JUNl) 

NRC=NRC-1 

SFRAC=SFRAC  -  FfiAC 

WR I TE( 6 . 36 ) I NCNUM , LUNL . LUNL 

WR I TE  <  7 , 36 ) I NCNUM , LUNL . LUN  L 

36  FORMAT ( //2X, 'ON  INCREMENT 14, '  ROLLER', 14, *  WAS  UNLOAOEO'/ 

2X, ‘ THIS  INCREMENT  NEGLECTED  AND  ROLLER 14, ' WILL  NOT*/ 
2X , ‘BE  CONSIDERED  IN  NEXT  INCREMENT') 

WRITE(7, 19) 

WRITE(7,12)( NODE( IP) ,GAR( IP ) ,  IP«1 , NR ) 

KEY  =  2 
GO  TO  109 
128  CONTINUE 

SFRAC=SFRAC  -  FRAC  +  ufrac 
00  130  IU=1 ,NR 
X I NC ( IU } =0 . 0 
FINC(IU)=0.0 
FT( IU)=0.0 

130  XT ( IU ) =0 . 0 
R  L  =  0 . 0 
RLINC=0.0 

DO  131  J= 1 , NR 
L=NODE( J) 

L  L  =  I RONO  < L) 

XINC(L)-UFRAC*BL(LL,1) 

XT(L)=X(L)  +  XINC(L) 

I F ( J.GE.NRC)GO  TO  132 
FINC(L)*-UFRAC*BL< LL,1)*R0LK 
FT(L)«F(L)  +  FINC(L) 

RL=RL  +  F ( L ) 

RIINC  =  RUNC  +  FINC(L) 

132  CONTINUE 

131  CONTINUE 
WRITE(6,37) I NCNUM, LUNL 
WR1 7  E ( 7 , 37 ) I NCNUM , LUNL 
GAPTEM=GAP( JUNL) 

DO  133  JJ=JUNL,NRC 
NODE ( J J ) =NODE ( J  J+1 ) 

133  GAP( JJ)=GAPI JJ+1 ) 

NODE ( NRC+ 1 )=NODE( JUNL) 

GAP(NRC+1 )=GAP(0UNL) 

NRC=NRC  -  2 

37  FORMAT ( //2X , 'ON  INCREMENT  NO. ',14,'  ROLLER  NO. ',14,'  WAS'/ 

.  2X, 'UNLOADED,  LOAD  FRACTION  TO  UNLOAD  ROLLER'/ 

.  2X, 'APPLIED  TO  DISPLACEMENTS  AND  LOADS  AND'/ 

2X, 'UNLOADED  ROLLER  REMOVED  FRON  MODEL') 

WRITE(7,19) 

W:TITE( 7 , 1 2)  (NODE(  IP)  ,GAP(  IP) ,  IP«1  , NR) 

121  CONTINUE 

DO  123  IA>1 ,NR 
X( I A)«XT ( 1A) 

123  F(IA)»FT(IA) 

WRITE(7 ,21 ) 
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PLATE 


1 . 
2. 

3. 

4. 

5. 
8. 

7. 

8. 
9. 

10. 
11  . 
12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21  . 
22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 
31  . 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 
41  . 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 
51  . 

52. 

53. 

54. 

55. 


C0MPILER(0IAG*3) 

SUBROUTINE  PLATE( A , B,D, XNU , NEL, I  PR I NT ) 

DIMENSION  EK( 12. 12,3), GK( 147,27), IN0DE(4) 

DIMENSION  B(3) ,HEAD(2) 

C 

C  GENERATE  ELEMENT  STIFFNESS  MARTIX  ANO 

C  GLOBAL  STIFFNESS  MATRIX 

C 
C 

C  READ  DATA  FOR  GENERATION  OF  ELEMENT 

C  STIFFNESS  MATRIX 

C 

ICK-0 

10  FORMAT ( ) 

DO  199  K*1 ,3 

CALL  PFE(A,B(K),D,XNU,EK(1 .1 ,K)) 

DO  1 98  1*1,12 

JB*  I  +  1 
DO  198  J* JB , 12 

198  EK(I,J,K)cEK(J,I.K) 

IF(IPRINT.EQ.O)CALL  WMTXt  EK(1,1,K),12, 12 .HEAD, 1 2 ) 

199  CONTINUE 

2000  FORMAT (/1X.6E15.6) 

C 

C  SUCCESSIVELY  ADD  THE  ELEMENT  STIFFNESS  MATRIX 

C  INTO  THE  GLOBAL  STIFFNESS  MATRIX 

c  oo  loop  on  the  element  numbers 

NUMR* 1 47 
NUMC=27 
I8W=NUMC 
LAD*0 

DO  100  IEL-1 , NEL 
C 

C  GENERATE  INDEX  FOR  SELECTION  OF  ELEMENT 

C  STIFFNESS  MATRIX 

C 

KELM*M0D(IEL.6) 

IF(KELM.LE.3)KSEL»KELM 
I F( KELM. EQ.4 )KSEL*3 
Ie(KELM.EQ.5)KSEL*2 
IF(KELM.E(J.0)KSEL=1 
C 

C  GENERATE  ELEMENT  NODE  NUMBERS  FROM  ELEMENT  NOS. 

C 

IN0DE{ 1 ) = I EL+LAD 
I  NODE ( 2 ) *  I  NODE ( 1 )  +  1 
INCDE( 3 ) = INODE( 2 )  +  7 
IN0DE(4)*IN0DE(  1  )  +  7 
L=MO0(IEL,6) 

I F ( ( L . EG . 0 ) . AND . ( I E  L . NE . 1 ) ) LAD*  L AO+1 
C 

C  START  DOUBLE  NESTED  DO  LOOP  ON  THE  NODES  OF 

C  THE  PLATE  ELEMENT  OR  ON  THE  3X3  PARTITIONS 

C  OF  THE  ELEMENT  STIFFNESS  MATRIX 

DO  101  1*1,4 


PUTE 


58. 

C 

ddS  ANO  IIB  ARE  THE  BEGINING  COLUMN  AND  ROM 

59. 

C 

indices  in  the  element  stiffness  matrix 

60. 

c 

dB  AND  18  ARE  THE  BEGINNING  COLUMN  ANO  ROM 

61  . 

c 

INDICES  IN  THE  GLOBAL  STIFFNESS  MATRIX 

62. 

c 

63. 

dB=(lNODE(d)  -  1)0 

64. 

I B*  (  I NODE  (  I  )  -  1)0 

65. 

ddB=( J-1 ) *3 

66. 

IIB-U-1  )*3 

67. 

c 

68. 

c 

DOUBLE  NESTED  DO  LOOP  IN  A  3X3  PARTITION  OF 

69. 

c 

THE  ELEMENT  STIFFNESS  MATRIX 

70. 

c 

71  . 

DO  1 03  11=1.3 

72. 

DO  104  JJ=1 ,3 

73. 

c 

74. 

c 

THE  FOLLOWING  FORTRAN  STATEMENT  WOULD  ADD  THE 

75. 

c 

THE  ELEMENT  MATRIX  TO  THE  GLOBAL  MATRIX  IF 

76. 

c 

WE  WERE  USING  FULL  SYMMETRIC  STORAGE 

77. 

c 

GK(  16+1 I  ,  JB+dJ)=GK(  IB+I I  tdB+Jd)+EK(  I IB+11  ,  ddB+dd.KSEL) 

78. 

c 

79. 

c 

HERE  WE  USE  BAND  SYMMETRIC  STORAGE  FOR  WHICH 

80. 

c 

THE  ROW  INDEX  RETAINS  ITS  TRUE  VALUE.  A  COLUMN 

81  . 

c 

IN  THE  BAND  SYMMETRIC  STORAGE  REPRESENTS  A 

82. 

c 

CO-DIAGONAL  ANO  ALONG  A  CO-DIAGCNAl  THE  DIFFERENCE 

83. 

c 

BETWEEN  THE  ROW  INDEX  AND  THE  COLUMN  INDEX  IS 

84. 

c 

CONSTANT  AND  IS  USED  AS  THE  COLUMN  INDEX  IN  THE 

85. 

c 

BAND  SYMMETRIC  STORAGE.  THE  COLUMN  INDEX  IS  ECUAl  TO 

86. 

c 

THE  BANDWIDTH  MINUS  THE  DIFFERENCE  IN  THE  ROW  AND 

87. 

c 

COLUMN  INDICES. 

88. 

c 

IR  ANO  IC  ARE  THE  ROW  AND  COLUMN  INDICES  OF  THE 

89. 

c 

GLOSAL  STIFFNESS  MATRIX  IN  BANO  SYMMETRIC  STORAGE 

90. 

c 

91  . 

IR=IB+I 1 

92. 

I C  =  J B  *  JJ 

93. 

IF( 1R.L7.IC)  GO  TO  104 

94. 

IRRrIR 

95. 

ICC=IB»  -  (IR  -  IC) 

96. 

181 

GK (  I  RR ,  1 CC )  =GK (  I RR  ,  ICC )  +  EK ( 1 1 8+ 1 1 ,  ddB+dd ,  KSEL ) 

97. 

104 

CONTINUE 

98. 

103 

CONTINUE 

99. 

102 

CONT INUE 

100. 

101 

CONTINUE 

101  . 

100 

CONTINUE 

102. 

WRITE' 3, 2001  )((GK( IP.dP) ,IP»1 ,147) ,dP=1  ,27) 

103. 

2C01 

FORMAT ( 2X , 6E20 . 1 2 ) 

104. 

RETURN 

105. 

END 
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1  . 
2. 

3. 

4. 

5. 

6. 

7. 

8. 
9. 

10. 
11 . 
12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21 . 
22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 
31  . 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 
41  . 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 
51  . 

52. 

53. 

54. 

55. 

56. 

57. 


C 

c 

c 

c 

c 


SUBROUTINE  PFE(A.B,D,XNU,EK) 

GENERATE  PLATE  FINITE  ELEMENT  STIFFNESS 
MATRIX 

A.BtPLATE  WIDTH  AND  LENGTH 
D.XNU.PLATE  STIFFNESS  AND  POISSION  RATIO 

DIMENSION  EK( 12,12) 

C=A/B 

C=C*C 

U*1 .0/C 

AS=A*A 

BS=B*B 

AB=A*B 

01=1.0  +5*XNU 
T3=1 .0/35.0 
T2=1 .0/25.0 

EX(1  ,1 )  =  156«T3*(C  +  U)  ♦  72»T2 
EK(4,4)  =  EK(1  .1  ) 

EK(7,7)*EK(1 .1) 

EK( ’0. 10 )=EK(1 ,1 ) 

EK(2,2)=(4*T3«C  +  52*T3*U  ♦  8*T2)»BS 
EK( 5 ,5) =EK(2 .2) 

EK(B,8)=EK(2,2) 

EK( l 1 . 11 )=EK(2,2) 

EK(3.3)=(52*T3*C  +  4*T3*U  +  8*T2)*AS 
EK(6.6)=EK(3.3) 

EK ( 9 , 9 ) =EK ( 3 . 3 ) 

EK(12,12)=EK(3,3) 

EK(2. 1 )=(22*T3*C  +  78*T3*U  ♦  6*T2rDl)*B 
EK(5,41»-EK(2,1) 

EK(B.7)*-EK(2.1  ) 

EM  1 1 , 1 0 )  sEK  (2 , 1 ) 

EK(3.1  )=-(78*T3.C  ♦  22*T3»U  ♦  6*T2*D1)*A 
EK(6,4)=EK(3.1) 

EK( 9 . 7 )=“EK( 3,1) 

EK( 12, 10)=-EK(3,1 ) 

EK(3.2)=-(11*T3«(C  +  U)  +  0. 5*T2*( 1 ,0+60*XNU) )»AB 
EK|6,5)=-EK(3.2) 

EK ( 9 , 8 ) =EK ( 3 . 2 ) 

EK ( 1 2 , 1 1 ) *-EK( 3 ,2) 

EK(4, 1 )=54*T3*C  -  156*T3*U  *•  72*T2 
EK{ 1 0, 7)=EK( 4, 1 ) 

EK(4.2)=( 13*T3*C  -  78*T3*U  -  6*T2)*B 
EK(5,1 )*-EK(4,2) 

EK(10,B)=-EK(4,2) 

EK( 1 1 ,7)>EK(4,2) 

EK(4,3)=(-27*T3*C  ♦  22»T3*U  +  6*T2*D1)*A 
EK ( 6 . 1 )=EK(4 ,3) 

EK(10.9)*-EK(4,3) 

EK ( 1 2 , 7 ) =-EK (4,3) 

EK( 5 , 3 ) * ( 13*0.5* T3*C  -  11*T3*U  -0 . 5*T2*D1 ) * AB 
EK(6,2)*-EK(5,3) 

EK( 1 1 , 9 ) »EK( 5 , 3 ) 

EK(12,8)«-EK(5.3) 

EK(5.2)«(-3*T3*C  ♦  26»T3*U  -  2*T2)«BS 
EK(11 ,8)«EK(5,2) 

EK(6.3)*(18*T3*C  -  4*T3*U  -  8*T2)*AS 
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58. 

EK( 1 2, 9 ) =  EK(  6 i 3) 

59. 

EK(7,1 )=-54*T3«(C 

♦  U)  +  72*T2 

60. 

EK(10.4)=EK(7,1 ) 

61  . 

EK( 8 , 1 ) * ( 1 3*T3*C  ♦ 

27*T3*U  -  6*T2)*B 

62. 

EK( 1 1 ,4 )*-EK(8, 1 ) 

63. 

EK(7.2)«- EK(8,1) 

64. 

EK( 1 3, 5 ) =EK ( 8, 1 ) 

65. 

EK(9. 1 ) = ( -27*T3»C 

-  13*T3*U  6*T2)*A 

66. 

EK (12,4)= EK (9,1 ) 

67. 

EK(7.3)=-EK(9,1 ) 

68. 

EK  ( 1 0 , 6 )  =— EK  (  9 , 1 ) 

69. 

EK ( 8 , 2 ) = ( 3*T3*C  + 

9*T3*U  +  2»T2)*BS 

70. 

E  K ( 1 1 , 5 ) =EK ( 8 , 2 ) 

71  . 

EK(9,3)=(9*T3»C  + 

3»T3*U  +  2*T2 ) *AS 

72. 

EK ( 1  2, 6 ) =EK( 9,3) 

73. 

EK(9,2)=(-13*0.5*T3*(C  +  U)  +0.5=T2)*AB 

74. 

EK(8,3)=+EK(9,2) 

75. 

EK( 12,5)=-EK(9,2) 

76. 

EK ( 1 1 ,6)=-EK(9,2) 

77. 

EK ( 1 0 , 1 )=-156«T3*C 

+  54*T3»U  -  72*T2 

78. 

EK (7,4)- EK (10,1) 

79. 

EK (11,1 )=(-22*T3«C 

+  27*T3»U  -  6«T2*D1)*B 

80. 

EK ( 1 0 , 2 )  =  EK ( 11,1) 

81  . 

EK(8,4)=-EK( 1 1 .1 ) 

82. 

EK(7,5)=-EK(11.1) 

83. 

EK ( 1 2 . 1  )  =  (-76*T3«C 

♦  1 3*T3»U  -  6*T2 ) * A 

84. 

EK(10,3)=-EK( 12,1 ) 

85. 

EK  ( 9 , 4 )  =  EK ( 1 2 , 1 ) 

86. 

EK(7,6)=-EK( 12. 1  ) 

87. 

EK( 11,2)- (-4*T2*C 

+  18*T3*U  -  8*T2)*BS 

88. 

EK(8.5)=EK(1 1 ,2) 

89. 

EK(12,2)=(-1 1*T3*C 

+  1  3*0 . 5*T3»U  -  0 . 5»T2*D1 ) *AB 

90. 

EK ( 1 1 , 3 ) =-EK (12,2) 

91  . 

EK ( 9 , 5 )  =  -EK ( 1 2 , 2 ) 

92. 

EK(6,6)=EK(12,2) 

93. 

EK ( 1  2 , 3  )  =  ( 26*T3«C 

-  3»T3*U  -  2*T2 ) *A5 

94. 

EK ( 9 , 6 ) =EK ( 1 2 , 3) 

95. 

DO  100  1=1,12 

96. 

DO  1 00  J=1  ,  X 

97. 

100  EK( I ,J)=EK( I ,J)'0/A3 

98. 

RETURN 

99. 

END 
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WRITEMATRIX 


1 . 

SUBROUTINE  WMT X( A , NR , NC .HEAD.NMAX) 

2. 

DIMENSION  A(NMAX.NMAX) ,HEAD(2) 

3. 

PRINT  202, HEAD 

4. 

K£  =  0 

5. 

KSET=NC/S 

6. 

KLEFT=r..CD(NC,8) 

7  . 

1F(KLEFT .NE.0)KSET*KSET+1 

6. 

DO  1  KT=1 ,KSET 

9. 

KB=KE+1 

1C. 

KE=kE+8 

1  1  . 

IF(KT.E0.KSET)KE=NC 

12. 

PRINT  200, (d,d=KB, KE) 

13. 

DO  1  1*1 ,NR 

14. 

1 

PRINT  201 ,1, (A(I,J).J«KB,KE) 

15. 

200 

FORMAT (//I  X, 10HROW  COL , 14, 7( 10X14) ) 

16. 

201 

FORMAT (1X,I4,8El4,6) 

17. 

202 

FORMAT ( ////I  X , 2A6 , 1 X ,6HMATRIX) 

18. 

RETURN 

19. 

END 

STSGAP2 


1  . 
2. 

3. 

4. 

5. 

6. 

7. 

8. 
9. 

to. 
11  . 
12. 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21  . 
22. 

23. 

24. 

25. 

26. 

27. 

28. 
29. 


SUBROUTINE  STGAP2(N0DE,G. L.GM.GSD) 

C  COMPUTE  THE  GAP  SIZES  BETWEEN  THE  ROLLERS 

C  AND  THE  PLATFORM  AS  A  SET  OF  RANDOMLY 

C  DISTRIBUTED  NUMBERS  OF  MEAN-GM  AND 

C  STANDARD  DEV  I  AT ION*GSD.  THERE  ARE  TO  BE 

C  L  OF  THESE  GAP  SIZES  IN  THE  ARRAY  G. 

C  THESE  GAP  SIZES  WILL  HAVE  A  GAMMA  DISTRIBUTION. 

C  THE  GAP  SIZES  ARE  ARRANGED  IN  ASCENDING  ORDER 

C  IN  ARRAY  G  AND  THE  CORRESPONDlN  ROLLER  NUMBERS 

C  ARE  IN  ARRAY  NODE. 

C  L-NUVBER  OF  ROLLERS 

C 

DIMENSION  G ( 1 ) ,NODE( 1 ) , WK ( 3) 

DOUBLE  PRECISION  DSEED 
LM3-  L 

CALL  ERTRAN(9, IR1 , IR2) 

'  DECODE ( 6 , 1 3 , I R2 ) JR 

13  FORMAT (16) 

DSEED*dR 
VARrGSD'GSD 
B* VAR/GM 
A=GM*GM/VAR 
WK( 1 )*0.0 

CALL  GGAMT(DSEED,A,B,LM3,WK,G) 

DO  1  1*1, LM3 

1  NODE ( I ) *  I 

CALL  A0RDER(G,N00E,LM3) 

RETURN 

END 
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AORDER 


x  xx  x  x 


1 . 

subroutine  AORDER( a, ia 

2. 

DIMENSION  A( N ) , I A( N ) 

3. 

LIMxN-1 

4. 

100 

I  NT  =  1 

5. 

DO  101  1=1 ,LIM 

6. 

IF(A(I+1 ).GE.A(I))  GO 

7. 

TEIUP*A( 1+1 ) 

8. 

I T  EM= I A ( 1  +  1 ) 

9. 

A( 1+1 ) =A( I ) 

10. 

I A ( 1  +  1  )  =  I A ( I ) 

1 1  . 

A; I jxTEMP 

12. 

I A  1 I )  =  I T  EM 

13. 

I  NT  =  I 

14. 

101 

CONTINUE 

15. 

I F ( 1NT.E0. 1 )  60  TO  102 

16. 

LIMxINT-1 

17. 

GU  TO  100 

18. 

102 

CONTINUE 

19. 

RETURN 

20. 

END 

•««•«  PLOAD  x=== 


1  . 
2. 

3. 

4. 

5. 

6. 

7. 

8. 
9. 

10. 
11 . 
12. 


SUBROUTINE  PL0AD(8L,P,NUMR0) 
DIMENSION  BL( 147,1) 

C 

C  P  IS  TAKEN  AS  THE  TOTAL  LOAD  SO 

C  P/NUMRO  IS  THE  L0A3  PER  ROLLER 

C 

r= P/NUMRO 
00  ICO  Ixl , NUMRO 
I P= I *3  -  2 
100  BL( IP.1 )=F 
RE1URN 
ENO 


(ROLLER)  STST/RUNST 


1  . 

•ASG.AX 

ECS'ORINT. 

2. 

OUSE  0 

.  .ECS'ORINT . 

3. 

•8RKPT 

PRINTS/O 

4. 

•PRT.S 

ECS*ROUER .  STST/RUNST 

5. 

OUSE  R ■ 

.  ,ECS*ROLLER. 

6. 

•ASG.T 

3. 

7. 

•ASG.AX 

ecs*output . 

8. 

•USE  7. 

,ECS*OUTPUT. 

9. 

•FOR , S 

R.MAIN2.M 

10. 

•MAP. IS 

A, 8 

11  . 

IN  M 

12. 

IN  R. 

13. 

LIB  DAO* IMSL . 

14. 

END 

15. 

•XQT  8 

16.  10.0,17.7,18.0,14.8,1630000.0,0.3,10000.0,150000.0 

17.  1,36.8 

18.  10.28,0.002,0.002 

19.  66RKPT  PRINT* 
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APPENDIX  D 


Integration  of  the  NASTRAN  Platform  Models  with 
the  Statistical  Model 

The  statistical  model  described  in  this  report  is  incomplete  in  that  it  is  missing  the 
description  of  the  platform-roller  interaction.  Should  data  be  obtained  to  resolve  this  knowledge 
gap,  it  would  be  desirable  to  integrate  the  NASTRAN  platform  structural  models  with  the 
program  for  carrying  out  the  solution  of  the  statistical  model  described  in  Appendix  C.  To 
carry  this  out  two  areas  of  activity  are  required:  modification  of  the  solution  algorithm  MAIN2 
given  in  Appendix  C,  and  obtaining  the  platform  stiffness  matrix  from  NASTRAN.  The  first 
of  these  areas  is  straightforward  programming  and  would  include  expanding  the  size  of  the 
arrays,  changing  the  input,  changes  to  account  for  the  rollers  that  are  offset  from  nodes  in 
the  NASTRAN  models  and  the  addition  of  an  out-of-core  equation  solver  to  handle  the  increased 
problem  size.  The  second  area  of  activity,  obtaining  the  platform  stiffness  matrices  from 
NASTRAN,  would  also  seem  straightforward  but  in  the  attempts  made  as  a  part  of  this  work 
it  turned  out  not  to  be.  The  NASTRAN  documentation  describes  auxiliary  subroutines  for 
the  output  of  data  such  as  stiffness  matrices  and  procedures  to  alter  the  rigid  formats  to  include 
these  subroutines.  However,  the  subroutine  that  appeared  to  provide  the  most  direct  means 
for  doing  this,  OUTPUT2,  did  not  work  properly  either  in  trying  to  write  information  on 
an  assigned  file  or  to  the  punch  file.  The  alternative  technique  used  was  to  insert  the  subroutine 
MATPRN  in  the  static  analysis  rigid  format  to  print  the  stiffness  matrix  as  part  of  the  printed 
output  file.  To  get  around  having  to  repunch  too  large  a  quantity  of  data  back  into  the 
machine,  we  break  pointed  the  print  file  to  a  permanently  assigned  file.  After  the  execution 
of  NASTRAN  this  permanelty  assigned  file  could  be  edited  to  obtain  the  stiffness  matrix  in 
the  format  dictated  by  the  subroutine  MATPRN.  This  matrix  could,  however,  be  read  in 
this  format  and  converted  to  storage  in  the  banded  symmetric  format  required  by  program 
MAIN2  for  solution  of  the  statistical  model.  It  is  the  platform  stiffness  matrix  without  the 
attachment  of  any  springs  simulating  rollers  that  is  desired,  so  it  is  necessary  to  modify  the 
generation  of  the  bulk  data  deck  so  the  scalar  spring  elements  are  not  included.  This  process 
has  been  completely  carried  out  for  the  8-foot  platform  and  has  been  read  into  a  modified 
version  of  MAIN2.  However,  debugging  of  the  modified  version  has  not  been  completed. 


52 


